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ABSTRACT 


This is a theoretical study to determine a practical design 
formula for the buckling load of simply supported rectangular plates 
under combined axial load and normal pressure. 


The theory and method used was based on Samuel Levy's 
integrated golution of the non-linear differential equations as de- 
rived by von Karman. The length-to-width ratio (a/b) of the plates 
considered was extended up to a value of 4:1. For each lengta-to- 
width ratio a set range of normal pressures was used in obtxining 
the different buckling loads. 


The study confirmed the conclusion arrived at by Levy; that 
is, the normal pressure increases the buckling load. It was further 
observed and concluded that for low values of a/b (less than 1.0) the 
buckling loads obtained in the manner described in Chapter III are 
impractical to attain. Yield failure will probably occur before the 
calculated critical load is reached. 


In general, however, the authors’ recommended buckling load 
formula for a steel plate is of the form, 


_ #°E rn? 
= sabia? * b 
where the value of K is scaled from Fig. XXXV, The basic limitation 
on the accuracy of this value of K is the approximated expression(s) of 
the deflection equation. The computed value of co, must be compared 
with the yield strength of the mate’ ial since the validity of the theory 
only holds true within the elastic region. 
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I. INTRODUCTION 


Background 


The discussion of buckling of plates is usually based upon a 
linear differential equation derived under the assumption that the 
deflections of the plates are small in comparison with their thickness. 
In particular, present design practice on hull plating under combined 
normal pressure and axial load is often based on approximate critical 
buckling load formulas. The normal load has generally been relegated 
to minor roles. The reasons purported to supported this were very well 
founded, both from the standpoints of theory and practice. Samuel 
Levy and others [2], [3} arrived at the conclusion that normal 
pressure increases the buckling load and thereby indicated that it 
would be conservative design to neglect the effect of normal pressure. 
The same conclusion was stated by Bleich {1]. Bleich further made 
mention of the fact that deflection of ship's hull plating usually does not 
exceed one-half the plate's thickness and therefore concluded that 
linearized theory could be applied. 


It is worthy to nzote at this point, however, that deflections 
exceeding half the plate's thickness are possible, especially in the 
case of simply supported plates. In those cases linear theory of 
plates no longer applies. As will be briefly explained in the sub- 
sequent paragraph, the problem becomes a non-linear stress problem 
since deflections of the order of magnitude of the plate thickness must 
be considered. 


Briefly stated, the nonlinearity of the differential equations has 
its origin in the fact that, in the case of large deflections, there is an 
interaction between the membrane stresses and the curvature of the 
plate. This interaction leads to non-linear terms in the eq'ations of 
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equilibrium of the plate elements. 


The complete differential equations of the problem were 
formulated by von Karman, who added those non-linear terms 
pertaining to the flexural rigidity of the plate. Timoshenko [6], 
Marguerre and Trefftz [8] derived the expressions for the strain 
energy of plates with large deflections. 


An attempt to develop the large deflection theory for rec- 
tangular piates under combined bending and longitudinal compression 
was made by Bengston {7}. His studies included plates with simply 
supported and clamped edges. However, Bengston introduced, in 
the course of his analysis, certain arbitrary assumptions which in 
part contradict each other and it is very doubtful whether the results 
of his analysis can be considered as entirely correct. 


In a series of papers Samuel Levy gave exact theoretical 
solutions for rectangular plates with large deflections. The 
theories developed include simnly supported and clamped plates. 


The problem of rectangular plates carrying longitudinal 
compression and normal! pressure is of prime importance in the 
design of the hull plating of ships. It is closely related to the 
question of buckling strength of plates. 


Statement of the Problem 


It is common knowledge that factors of safety are actually 
factors of ignorance which in some cases have remained unchanged 
for a period of time for no obvious reason. Factors of safety may 
be used to account for one or more of the following: (a) material 
imperfections, (b) faulty construction practices, (c) material de- 
gradations due to corrosion, and (d) lack of ability to make rigorous 
calculations. The conservativeness of design as previously mentioned 
prompted the authors to ask in what way would conservative design 
affect the factor of safety; that is, does conservative design actually 
compound the factor of safety? If, for instance the safety factor is 
increased by only a few tenths of a percent, then one is justified in 
using linear theory. However, if the increase becomes considerable, 
the design may become uneconomical. 
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With the salient points stated in the previous paragraph in 
mind and also with the conclusion of Levy on the effect of normal 
pressure on the buckling load of plates, the authors came to the 
conclusion that there is a need for further investigations. First, 
these investigations afford the possibility of checking the accuracy 
of simpler approximate methods and formulas. Secondly, and 
actually the ultimate goal of this thesis, these investigations would 
hopefully lead to a formulation of a simple design formula for plates 
subjected to the simultaneous action of normal and edge loads. 


A literature survey disclosed there is no simple formula which 
a designer can use with great facility. Fortunately, general solutions 
to the differential equations are available in the form of Fourier series. 
Though the analysis of these solutions is highly involved and the numeri- 
cal work for obtaining special solutions is very laborious, a method 
(Chapter III) has been devised which lends itself to computer programming. 


This thesis is concerned mainly with plates simply supported 
on all edges and subjected to combined loadings of uniform normal 
pressure and axial load in one direction. It is further hoped that it 
will serve as a forerunner of later investigations based on different 


boundary conditions. 
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Nomenclature (See Fig. I): 


a plate lemeth in the x-direction. 
b plate width in the y-direction. 
a alb 
h plate thickness. 
Ww plate deflection. 
ony, coordinate axes with origin at corner of plate. 
E Young's modulus. 
Ll Poisson's ratio. (Note: Tabulated results and figures are 
for steel with uw = 0. 300) 
uniform normal pressure on plate. 
average compressive strain at edges y = 0, b. 
axial load on plate. 
stress function. 
= Eh? /12 (1 - en Flexural rigidity of the plate. 
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Fig. I Schematic diagram of plate under axial load 
and normal pressure. 





Fundamental Equations 


The fundamental partial differential equations governing the 
deformations of flat plates have been derived by von Karman. They 
are given by Samuel Levy [2] and Timoshenko | 6}. The exact math- 
ematical analysis to dete: mine the buckling strength of a simply sup- 
ported flat plate under combined edge compression and normal loading 
involves the integration of these “rez yp Pies, 


4 4 44 be 
og + 2 ys ee yy | () 
x CY 


a4w » 9 2_W a *w ,h/2°F aew 9°F aw 2 oF ows (2) 
) 
py? ox? ax” ay? OY Oxy. 


The boundary conditions for a simply supported rectangular 
plate to be satisfied by equations (1) and (2) are for deflection, w, 
and edge bending moments per unit length to be zero at the edges 


of the plate; viz., ' 


‘ ac 
moa - DCS ay SH) - 0 at x= 0, x= a (3) 
Ox by 
y dew aw 
mye - DoS N+ yl) 20 at y20, yeb (4) 
J * Oy ox 


The complete solution is more fully explained in reference { 2] 
and, without loss of gene: ality, the deflection equation is here only 
reproduced in the —— Fourier series, 


marx nv 
Sy Wn <. sin - sin 17 (5) 
m=1 nx! 


where the undetermined constants, w. ., must satisfy the relation(s) 


expressed by equation (9) of the same reference. 


The average compressive stain. e, at the edges y = 0, b is 
computed from equation (11) of [ 2] as: 
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7 DOF 


In this thesis, the deflection equation for large values of a/b 
(2 >> 1) is approximated by the following expression, 


3 7% 7y Sax Ty 
Ww Way sin =~ sin - + ™., i sin . sin 5 
. OFX Ty | (mk a. Wy 
+ Ws 4 sin sin - rs Wo 4 gin — sin= , (7) 


and for a/b - values in the vicinity of 1 and less than 1, the approximated 
deflection equation is, 


"= Ww sin 1 gin S2¥ 


w7X oO. Ty 
ae + 
sin = sin 5 Ww 5 


1,1 1,3 


. on xX wy oa x say 
+ Ws , sin — sin 5 + We g Sin . sin 5 


WX oO. OMY, om x Ty 
+ w sin _ sin = +W sin - sin _ (3) 


1,5 5, 1 


There are no special reasons for using two deflection equations 
for different ranges of a values, except that equation (7) restricts the 
deflected shape of the plate to one sine wave across the width and a 
combination of four sine waves along its length. The errors involved 
by this approximation are expected to be less than five percent { 2}, 
{3}. For this reason, equation (7) could be reduced to a fewer 
number of terms as @ is reduced. This is evidenced by the deflection 
equation used fo: a = 3.0 in refe:xence [2]. More specifically, the 
contribution of the higher-ordered deflection coefficients in the x- 
direction becomes less significant as a is reduced from ane >> 1. 


However, as @ is further reduced the assumption of one sine 
wave across the width of the plate becomes incorrect since the 
higher-ordered deflection coefficients in the y-direction become 
increasingly significant and cannot be neglected unless the desired 
degree of accuracy is sacrificed. Thus, for suche values, equation 
(8) is used. As before, it is believed that the errors incurred by 
uging a finite number of deflection coefficients would give results 
within the desired accuracy (see Table 13 of [ 2]}. 
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The determination of the value of @ applicable to equation (7) 
or (8) and the more general case of expanding the work of Levy et al. 
on simply supported rectangular plates within the practical values of 
a/b are problems examined in this thesis. 





Il, PROCEDURE 


With some additional steps and explanatory comments incorporated, 
the procedure followed in solving for the deflection equations was essentially 
the same as that outlined in reference [3]. The steps were ae follows: 


1. The family of four or six simultaneous cubic equations 
corresponding to the same number of unknown deflection coefficients 
was first obtained. This was done by solving for the coefficients a. q? 
defined by equation (8) of reference [2], and substituting them into 
equation (9) of the same reference. The results are shown in Appendix A. 
3 


2. Each of the resulting equations was divided by h. This step 


nondimensionalized the family of equations. 


3. Velues of w, ,/h, wy 3/h, We ,/h, etc. were estimated 
corresponding to rare, aes of Pb/Eh? ‘and pb4/ Eh* -- the non- 
dimensional axial load and normal pressure, respectively. This was 
the most delicate step since at higher values of Pb/ Eh3 (normal pressure 
held constant) more than one solution was found possible (see { 3] also). 
In other words, the estimated w's, when not properly chosen, could lead 
to solutions other than those corresponding to a continuous chenge in 
buckled form from zero axial load to the buckling load. 


It was therefore decided that a step-by-step method of estimating 
the w's be employed. In conjunction with the subsequent steps 4, 5, and 
6, this method proceeds as follows: 


a). Corresponding to three low values of Pb/ Eh°, including 
Pb/Eh® =x 0, three sets of solutions were obtained using a reasonable 
set of w-estimates. These values of Pb/Eh* were hoped to lead 
to unique solutions. For these estimated w's, reference | 2] or 


{ 3] was used as a guide. 
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b). From (a) above, succeeding estimates of w's as 
axial load is increased, were based on the continued trend 
of each w-ve+Pb/Eh® curve; viz., the slope and the rate of 
change of slope computed from the last three correct points 
were employed in an abridged Taylor series of the form, 


i 
’ h ” h 
w(x, +h) = w(x,) + w (xp) itt w (x5) T° 


where, wx, + h) = next estimate of the particular w. 
w(x,) = last computed value of w. 


w'(x)) = slope obtained from the last two computed 
points. 


w’ (xp) rate of change of slope based on the last 


three computed values of w. 
4, The resulting cubic equations shown in Appendix A were 
linearized by expanding the right-hand side of each equation in a Taylor 


series in the neighborhood of the estimated values of Ww, 1/ h, w 3/ h, 


1, 
wa ,/h, ete., omitting higher ordered terms. 

5. Crout's method [5] was used in solving for the difference 
between the estimated w's and their improved values. 


6. Step 5 was repeated until such time tht the calculated 
error sum was less than atest constant. The error sum was defined 
by the authors ag the measure of closeness of successive approximations. 
In equation form, it is as follows; 

N 
Error = . x; 
rae, 
i= | 
where, x; represents the absolute value of the difference between an 
estimate w and its improved value. 


The test constant used for equation (7) was 0.0001 and for 
equation (8), 0.000001. Why two different values of the test constant 
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were used would become evident from a study of the results in 
Chapter IV. 


7. The ave-age strain corresponding to each axial load was 


computed from equation (6). 


8. Inthe vicinity of the buckling load where the strain started 
to change more rapidly than the axial x cg the Why /h was made the 
independent variable in place of Pp /Eh?. The other deflection coef- 
ficients were computed in a similar fashion as before. 


Values of the buckling load were obtained for values of pb* / Eh* 
equal to 2. 50, 7.50, 12.50, 18.00, 24.50, and 30.00, and at values of 
a/b up to 4.00. These were plotted versus a/b and analyzed for a 
possible comparison with Bryan's classical solution. 


The whole procedure(s) had been programmed with the use of 
the IBM-7094. Details not oth>rwise covered here may be found in 
the programs shown in Aprendix B. 
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IV. PESULTS 


Tables Ito XII are sample results of the computer programs. 
Tables I to VI were obtained using equation (8) for an a/b = 1.00 and, 
the rest of the tables were for an a/b = 4.00 using equation (7). 


Figures II to XXXIII are the graphical results of the same 
programs. The curves drawn are the stress-strain curve analogy 
for plates with the non-dimensional axial load, Pb/ Eh, as ordinate 
and the non-dimensional average strain, eb” / hn’. as abscissa. (Note: 
Not all results for the chosen a/b-values used were plotted. ) 


Tabulated values of the critical loads, rounded to the nearest 
hundredth, are shown in Tables XII and XIV. 


Figure XXXIV is a plot of the critical Pb/ Eh? versus a/b. The 
curves between zero a/b and the smallest value of a/b used have been 
extrapolated. The curves have been faired at points where the results 
of both equations (7) and (8) either become tangent or intersect. Re- 
finements of the curves in the vicinity of the cusps were made by ob- 
taining more results at the a/b-values concerned. These results are 
shown in Tables ATI and AIV. 


= ¢ a = —— _— , 

: . : ' 
6 

7 
) . . a } 
 -« — | | 

0 
| > 

« " | 

| : ° * A 

2 ron y = 

we « “i - 


Moria oyrreve lanoiue ied 
| aren ae 


— . 
asi ; ha, , : , . ~ 
| | | * 








12 


Table I - Values of deflection coefficients for various values of axial compressive load in the 
x-direction, P, for simply supported rectangular plate, a/b = 1.00, u = 0.300. 
Normal pressure, p = 2. 50 Ent /b?, 
Pb “1,4 “1,3 "3,1 "3,3 "1,5 bd eb 
Fhe h h h h h h ho 
0. 0.113109 0.001522 0.001522 0.000156 0.000134 0.000134 ‘0.015 
0.500 0.130955 0.001535 0.001606 0.000159 0.000135 0.000137 0.521 
1.000 0. 155300 0.001552 0.001704 0.000162 0.000135 0.000140 1.029 
1. 500 0.190218 0.001580 0.001826 0.000166 0.000135 0.006143 1.544 
2.000 0.243439 0.001636 0.001995 0.000170 0.000135 0.000146 2.073 
2. 300 0.329375 0.001781 0.002298 0.000177 0.000135 0.000150 2.633 
3.000 0.466786 0.002213 0.003028 0.000195 0.900136 0.000154 3.268 
3. 500 0.651813 0.003327 0.004820 0.000246 0.000137 0.000161. 4.024 
4,000 0.848735 0.005355 0.008196 0.000384 0.000143 0.000176 4.889 
4,500 1.034310 0.008183 0.013194 0.000667 0.900154 6.000205. 5. 821 
9. 000 1.204318 0.011610 0.019678 0.001151 0.000172 0.900252 6. 793 
5. 500 1. 360460 0.015481 Q., 027524 0.001886 0.000196 0.000320 7.792 
6,000 1.505171 0.019690 0,036648 0.002916 0.000227 0.000412 8.810 
6. 500 1. 640526 0.024166 0.046998 0.004283 0.000264 0.000530 9.845 
7,000 1. 768140 0.028864 0.058539 0.006030 0.000308 0.000674 10. 896 
7. 500 1.889255 0.033752 0.071252 0.008198 0.000359 0.000846 11.961 
8.000 2.004827 0.038810 0.085125 0.010828 0.006420 0.001046 13.042 
8. 500 2.115604 0.044023 0.100152 9.013961 0.000492 0.001273 14,137 
9.000 2.222177 0.049379 0.116321 0.017636 0.000580 0.001528 15. 248 
9. Sad 2.325014 0.954867 0.133058 0.021890 0.000688 0.001307 16.376 
10. 000 2.424491 0.960477 0.182120 0.026757 0.000821 0.002110 17, 521 
10. 500 2.520910 0.066194 0.171741 0.032266 0.000986 0.002434 18. 684 
11.000 2.614510 0.071999 0, 592479 0.038442 0.001190 0.002775 19. 867 
11, 500 2.705480 0.077869 0.214327 0.045301 0.001442 0.603128 21.070 
12,000 2.793967 0.083776 0.237259 0.052850 0.001751 0.003489 22,295 
12, 500 2.880085 0.089685 0. 261241 0.061089 0.002128 0.003851 23.543 
13.000 2.963914 0.095555 0: 286231 0.070006 0.002583 0.004206 24.813 
13. 500 3.045514 0.101342 0.312181 0,079581 0.003127 0.004546 26.108 
14, 000 3.124924 0.106995 0. 339036 0.089786 0.003771 0.004861 27.427 
14, 500 3.202167 0.112462 0. 366738 0. 100586 0.004525 0.005139 28.772 
15. 000 3.277257 0.117688 G. 395225 0.111941 0.005399 0.005370 30,142 
15. 500 3.350196 0.122617 0. 424455 0.123810 0.006404 0.005540 31.537 
16,000 3. 420982 0.127192 0.454368 9.136152 0.007549 0.005637 32.957 
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Table I (cont.) - Values of deflection coefficients for various values of axial compressive load in 
the x-direction, P, for simply supported rectangular plate, a/b = 1.00, u = 0.300. 


Normal pressure, p = 2.50 Eh*/b*. 











Pb af “1,3 = ¥I “3,3 "os eb? 
Fh? h h h h h h ne 
16.500 3.489605 0.131359 0.484926 0.148926 0.008843 0.005646 34.402 
17.000 3.556047 0.135060 0.516103 0.162099 0.010296 0.005552 35,873 
17.500 3.620286 0.138240 0.547881 0.175642 0.011916 0.005341 37.369 
18.000 3.682288 0.140841 0.580262 0.189534 0.013714 0.004996 38.891 
18.500 3.742007 0.142805 0.613260 0.203761 0.015701 0.004502 40.438 
19.000 3.799381 0.144069 0.646908 0.218318 0.017886 0.003841 42.011 
19.500 3.854328 0.144565 0.681258 0.233211 0.020285 0.002996 43.611 
20.000 3.906737 0.144218 0.716381 0.248451 0.022910 0.001948 45.239 
20.500 3.956468 0.142943 0.752370 0.264062 0.025780 0.000677 46.897 
21.000 4.003334 0.140647 0.789343 0.280075 0.028913 -0.000837 48.586 
21.500 4.047101 0.137221 0.827445 0.296530 0.032333 -0.002613 50.310 
22.000 4.087464 0.132543 0.866846 0.313473 0.036063 -0.004672 52.070 
22.500 4.124046 0.126482 0.997749 0.330955 0.040133 -0.007032 53.871 
23.000 4.156369 0.118897 0.950386 0.349025 0.044572 -0.009707 55.717 
23.500 4.183853 0.109654  v.995010 0.367720 0.049409 -0.012703 57.612 
24.000 4.205803 0.098647 1.041483 0.387049 0.054667 -0.016014 59.562 
24.500 4.221416 0.085840 1.091255 0.406970 0.060358 -0.019621 61.571 
25.000 4.229790 0.071321 1.143337 0.427361 0.066472 -0.023485 63.643 
25.500 4.229896 0.055357 1.198303 0.447995 0.072967 -0.027555 65.779 
26.015 4.219896 0.037904 1.258260 0.469164 0.079986 ~0.031938 68.049 
26.282 4.209896 6.028733 1.290778 6.479911 0.083709 -0.034291 69.250 
26.479 4.199896 0.021994 1.315644 0.487706 0.089491 -0.036091 70.151 
26.638 4.189895 0.016648 1.336290 0.493858 0.088743 -0.037593 70.885 
26.771 4.179895 0.012261 1.354111 0.498900 0.090636 -0.038904 71.506 
26.885 4.169895 0.008598 1.369846 0.503121 0.092259 -0.040080 72.043 
26.984 4.159895 0.005512 1.383944 0.506695 0.093668 -0.041154 72.513 
27.070 4.149895 0.002903 1.396704 0.509740 0.094900 -0.042148 72.927 
27.145 4.139895 0.000699 1.408341 0.512342 0.095983 -0.043078 73,295 
27.211 4.129895 -0.001158 1.419012 0.514563 0.096937 -0.043956 73.622 
27.270 4.119895 ~-0.002711 1.428842 0.516454 0.097777 -0.044790 73.914 
27.320 4.109895 ~-0.003996 1.437929 0.518054 0.098517 -0.045587 74.173 
27.365 4.099895 -0.005043 1.446354 0.519395 0.099167 -0.046352 174.404 
27.403 4.089895 -0.005877 1.454183 0.520505 0.099737 -0.047089 74.609 
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Table I(cont.) - Values of deflection coefficients for various values of axial compressive load in 
the x-direction, P, for simply supported rectangular plate, a/b = 1.00, u = 0. 300. 


Normal pressure, p = 2.50 B® fo*. 


Ph Wy “13 3,1 ¥3,3 ~s, “es eb 
Eh? - we h h he - wo “he 
27.437 4.079895 -0.006518 1.461472 0.521407 0.100238 -0.047802 74.790 
27.465 4.069895 -0.006986 1.468271 0.522120 0.100663 -0.048493 74.949 
27.489 4.059895 -0.007297 1.474621 0.522662 0.101031 -0.049164 75.088 
27.509 4.049895 -0.007465 1.480560 0.523047 0.101344 -0.045818 75.209 
27.524 4.039895 -0.007503 1.486119 0.523290 0.101605 -0.050457 75.313 
27.537 4.029895 -0,007422 1.491328  0.523401° 0.101819 -0.051080 75.402 
27,546: 4.019895 -0.007231 1.496211 0.523890 0.101989 -0.051691 75.475 
27.551 4.009895 -0.006940 1.500794 0.523268 0.102118 -0.052288 75,535 
27,554 3.999895 -0.006557 1.505095 0.523043 0.102210 -0.052874 75. 582 
27.555 3.989895 -0.006090 1.509135 0.522723 0.102266 -0.053449 75.617 
27.552 3.979894 -0,005545 i,512931 0.5223:4 0.102291 -0.054014 75.641 





Table II - 


Q. 

0. 500 
1, 000 
1, 500 
2.000 
2. 500 
3, 000 
3. 500 
4.000 
4.500 
2. 000 
9. 500 
6.000 
6. 500 
7. 000 
7.900 
8, 000 
8, 500 
9, 000 
9, 500 
10. 000 
10, 500 
11. 000 
11, 500 
12,000 
12. 500 
13.000 
13. 500 
14, 000 
14, 500 
15.000 
15. 590 
16. 000 


Values of deflection coefficients for various values of axial compressive load in the 


x-direction, P, for simply supported rectangular plate, a/b = 1.00, uw = 0.300. 


Normal pressure, p=7. 50 En! /p*. 


1,1 
h 


0. 329069 
0.375274 
0.433827 
0. 508237 
0.601481 
Q.714013 
0.842155 
0.979252 
1.118822 
1, 256523 
1.390149 
1, 518869 
1.642574 
1. 761487 
1.875953 
1.986346 
2.093025 
2.196316 
2.296504 
2. 393836 
2.488521 
2.580733 
2.670615 
2.758281 
2.843819 
2.927293 
3.008751 
3.088220 
3. 165716 
3.241245 
3.314803 
3. 386382 
3. 455968 


dV 
h 


0.004721 
0, 004842 
0.005031 
0.005343 
0, 005870 
0.006741 
0, 008089 
0.009994 
0.012457 
0,015424 
0.018820 
0.022574 
0.026629 
0. 030937 
0.035466 
0.040189 
0.045088 
0.050148 
0.055356 
0, 060701 
0.066170 
0.071745 
0.077408 
0.083135 
0.088897 
0.094660 
0.100383 
0, 106024 
0.111533 
0.116862 
0.121956 
0.126763 
0.131229 





0.004721 
0.005064 
0.005516 
0.006154 
0.007116 
0, 008610 
0.010891 
0.014177 
0.018599 
0. 024205 
0.030995 
0. 038954 
0. 048062 
0. 058305 
0.069672 
0.032157 
0.095759 
0.110476 
0. 126310 
0. 143262 
0. 161328 
0. 180504 
0. 206779 
0, 222136 
0.244550 
0. 267990 
0.292416 
0.317781 
0.344034 
0. 371120 
0. 398934 
0.427576 
0.456848 


0.000485 
0.000499 
0.000518 
0.000545 
0.000587 
0.000658 
0.000781 
0, 000988 
0.001321 
0.904821 
0,962531 
0.002492 
0.004746 
0.006331 
0.008287 
0.010655 
0,013473 
0.016780 
0.020613 
0.025007 
0.029994 
0, 035602 
0.041851 
0.048757 
0.056327 
0.064559 
0.073442 
0.082957 
0.093078 
0. 103773 
0, 115003 
0.126731 
0.138918 


~~. 5 


h 





0.000403 
0.000404 
0.000404 
0.000405 
0.000406 
0.000408 
0.000411 
0.000418 
0.000429 
6, 000446 
0, 000467 
0,000494 
0.000527 
0, 000566 
0, 000612 
0.000667 
0.000732 
0, 000809 
0.000904 
0.001019 
0,001160 
0.001335 
0.001550 
0.001813 
0.002133 
0.002529 
0.002985 
0.003536 
0.004186 
0.004942 
0. 005816 
0.006817 


0.007954 


0.000403 
0.000412 
0.000421 
0.095430 
0.000441 
0.000454 
0.000472 
0.000496 
0.000531 
0.000580 
0.000646 
0.000731 
0. 000839 
0. 000970 
0.001126 
0.001308 
0.001515 
0.001749 
0.002007 
0.002290 
0.002594 
0.002919 
0.003260 
0.003613 
0.003973 
0.004334 
0.004688 
0.005028 
0.005344 
0.005625 
0.005860 


0. 006037. 


0.006143 


eb" 
- 


0,133 
0.674 
1,232 
1.819 
2.446 
3.129 
3.876 
4,685 
9. 948 
6, 454 
7, 395 
8. 363 
9,355 
10, 367 
11.397 
12,445 
13.510 
14,593 
15. 692 
16. 809 
17,944 
19.099 
20.273 
21. 469 
22. 686 
23.927 
25.190 
26.478 
27.790 
29.127 
30, 489 
31.876 
33. 289 
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Table II (cont.) - Values of deflection coefficients for various values of axial compressive load in 


Pb 


tm 


Fh? 


16. 500 
17.000 
17. 500 
18. 000 
18. 500 
19. 009 
19. 500 
20.000 
20. 500 
21.000 


21. 500 
22.000 


22.500 
23.000 
23. 500 
24.000 
24. 500 
25.000 
25. 500 
25.550 
25. 600 
26. 343 
26. 641 
26.859 
27,035 
27.183 
27, 309 
27,419 
27.514 
27,599 
27. 673 
27. 738 
27, 796 


the x-direction, P, for simply supported rectangular plate, a/b = 1.00, u = 0.300. 
Normal pressure, p = 7.50 Eh‘ /b : 


Ya 


1, 
h 


3.523541 
3. 589080 
3.652556 
3.713931 
3. 776869_ 
3. 850180 
3.884915 
3.937264 
3.587097 
4,03454b 
4.078508 


4.119609 
4.157217 
4.190920 
4.220217 
4.244515 
4.263144 
4.275380 
4.280478 
4.280566 
4. 280573 
4.270575 
4.260573 
4.250573 
4.240574 
4.230573 
4, 220573 
4.210573 
4.200573 
4.190573 
4.180573 
4.170573 


4.160573 


7.3 


wa 

0.135299 
0.138921 
0.142039 
0.144599 
0.146544 
0.147813 
0. 148342 
0. 148057 
0.146879 
0.144717 


0.141468 
0.137016 


0.131232 
0.123980 
0.115123 
0.104545 
0.092182 
0.078069 
0.062397 
0.060758 
0.059106 
0.033610 
0.023209 
0.015668 
0.009719 
0.004853 
0.000801 


-0. 002606 
-0.005481 
-~0. 007906 
-0.009946 
-0.011650 
-0. 013058 


W341 


A 





0. 486762 
0.517293 
0. 548417 
0. 580139 
0.612469 
0. 645434 
0, 679081 
0.713474 
0. 748695 
0.734851 


0. 822070 
0. 860509 


0.900349 
0.941796 
0.985077 
1.030425 
1, 078057 
1.128141 
1.180770 
1.186174 
1.191604 
1.275424 
1.310761 
1. 337450 
1. 359503 
1. 378486 
1, 395211 
1.410168 
1, 423684 
1. 435988 
1, 447253 
1. 457613 
1.467173 


4 


“3,3 


a 


h 


0.151527 
0.164526 
0.177889 
0.191594 
G. 205629 
0.219990 
0.234680 
0.249711 
0.265104 
G, 289889 
0, 287103 
0. 313792 
0. 331006 
0.348794 
0, 367198 
0. 386239 
0, 405896 
0.426081 
0.446614 
0.448676 
0.450739 
0, 481195 
0, 493098 
0.501628 
0, 508332 
0.513818 
0.518406 
0.522289 
0.525598 
0. 528424 
0.530839 
0.532896 
0.534639 


71,5 


h 


0. 009236 
0.010672 
0.012271 
0.614042 
0.015995 
0.013142 
0.020494 
0.023066 
0.025874 
0.028936 
0.032274 
0.035911 
0.039873 
0.044191 
0. 048890 
0. 053997 
0.059526 
0. 065474 
0.071809 
0.072461 
0.073117 
0.083174 
0.087319 
0. 090380 
0.092848 
0.094917 
0.096689 
0.098225 
0.099567 
0.100745 
0.101782 
0.102695 
0. 103498 


Ws 1 


nin ee 


oF 


0.006165 
0. 006088 
0.005397 
0.005579 
0.005116 
0.004493 
0. 003692 
0.002695 
0.001484 
0, 006040 


-0. 061656 
-0. 003627 
-0. 005885 
-0, 008447 
-0.011318 
-0. 014494 
-0. 017955 
-0. 021661 
-0, 025553 
-0, 025950 
-0,. 026348 
-0.032375 
-0. 034850 
-0. 036706 
-0, 038241 
-0.039573 
-0. 040763 
-0, 041847 
-0. 042849 
-0.043785 
-0.044669 
-0,045509 
-0.046313 


34. 726 
36.188 
37. 676 
39.188 
40.726 
42. 289 
43.679 
45. 497 
47,143 
48.6821 
50. 531 
52.277 
94. 062 
95. 891 
97. 767 
99.695 
61. 680 
63.725 
65. 831 
66.045 
66.259 
69.518 
70. 860 
71.856 
72.667 
73.353 
73.946 
74.465 
74,924 
75. 332 
75. 694 
76,018 
76. 303 
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Table If (cont.) - Values of deflection coefficients for various values of axial compressive load in 


27. 846 
27.890 
27.928 
27.960 
27. 988 
28.011 
28.030 
28,045 
28.057 
28.065 
28.070 
28.072 
28.071 


ae 


b 


4.150573 
4.140573 
4.130573 
4.120573 
4.110573 
4.100573 
4.090573 
4.080572 
4.070572 
4.060572 
4.050572 
4.040572 
4.030572 


1.3 


h 


-0.014204 
-0.151116 
-0.015817 
-0. 016330 
-0.016670 
-0.016855 
-0,016899 
-U. 016813 
-0. 016609 
-0.016298 
~0. 015887 
-0.015385 
-0,.014799 


W301 
h 


1.+76021 
1. 484228 
1. 491858 
1.498961 
1, 505584 
1. 511767 
1.517544 
1. 522947 
1, 528004 
1. 532740 
1.537178 
1. 541339 
1. 545240 


Normal pressure, p = 7. 50 Fh‘ fp‘. 


0.536103 
0.537317 
0. 538308 
0.539097 
0.539701 
0.540138 
0. 540422 
0.540566 
0. 540580 
0. 540476 
0. 540262 
0. 539947 
0. 539538 


0.104203 
0.104821 
0. 105358 
0.105823 
0.106222 
0. 106560 
0.106842 
Q. 107073 
@.107257 
0.107397 
0.107497 
0.107559 
0.107586 


-0.047085 
-0. 047830 
“0. 0498557 
-0. 049252 
-0. 049934 
-0. 050600 
-0. 051251 
-0. 051888 
-0.052513 
-0, 053127 
-0. 053730 
-0.054322 
~0, 054905 


the x-direction, P, for simply supported rectangular plate, a/b = 1.00, = 0.300. 


76. 563 
76.791 
76.993 
T7171 
77. 328 
77. 464 
77. 582 
77. 683 
77. 767 
77, 837 
77, 892 
77.935 
77.965 
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= Table Ht ‘ Values of deflection coefficients for various values of axial compressive load in 


the x-direction, P, for simply supported rectangular plate, a/b=1.00, :*0. 300. 


Normal pressure, p# 12.50 Eh*/b A 


4 





Pb * 4 As W3 1 W3 3 Wiis "5.1 eb” 
> —-— «= . ge ee 
0. 0.521574 0.008252 0.008252 0.000856 0.000672 0.000672 0, 336 
.500 0.583586 0.008596 0.008983 0.000897 0.000673 0.000687 0,921 
1,000 0.656741 0.009091 0,009947  0,000951 0.000674 0.000702 1,533 
1.500 0.741808 0.009806 0.011249  0,001029 0.000676 0.000720 2.180 
2,000 0.838433 0.010816 0.013027 0.001143 0.000679 0.000740 2. 869 
2.500 0,944901 0.012197 0.015433 0.001312 0.000684 0.000764 3. 604 
3.000 1.058458  0,013999 0.018613 0.001562 0.000691 0.000795 4. 386 
3.500 1.176059 0.016239 0.022681  0,001923  0,00070i 0.000835 5.212 
4.000 1.295044  0,013902 0.027715 0.002429  0.00071G 0.000886 6. 078 
4.500 1.413472 0.021955 0.033760  0,003116 0.000733 0.000951 6,978 
5.000 1.530103 0,025354 0.040840 0.004021 6.000760 0, 001032 7,907 
28,582 4.108265 -0,025822 1.565221 0.557478 0.112708 - 0.053915 80, 164 
28.589 4.098265 -0,025366 1.569745 0.557257 0.112808 -0,054531 80.226 
28.592 4.088265 -0.024813 1.573976 0.556930 0.112869 -0.055189 80.274 
28.592 4.078265 -0,024172 1.577937 0.556504 0.112893 -0.055739 80. 308 
28.590 4.068265 -0.023451 1.581643 0.555988 0.112882 -0.056330 80.321 





Table IV 


tS 
on 
© 
© 


29.156 
29.164 
29.169 
29.170 
29.169 


- Values of deflection coefficients for various values of axial compressive load in the 


x-direction, P, for simply supported rectangular plate, a/b = 1.00, u 


Normal pressure, p= 18.00 Eh*/b?. 


0. 704259 
0.773978 
0.852085 
0.938231 
1.031473 
1. 130381 
1, 233292 
1.338575 
1. 444328 
1, 550963 
1. 656205 


4.159187 
4.149187 
4.139187 
4.129187 
4.119187 


Wi 3 
h 





0.012496 
0.013128 
0.013965 
0.015057 
0.016453 
0.018189 
0.020282 
0.022732 
0.025526 
0.028639 
0.032044 


-0, 036514 
-0.036003 
-0. 035390 
-0. 034684 
-0.033893 


ws; 
h 


0.012496 
0.013706 
0.015241 
0.017200 
0.019688 
0.022809 
0, 026658 
0.031311 
0.036825 
0,043241 
0.050589 


1.600875 
1.605478 
1, 609775 
1.613785 
1.617529 


“3,5 
“h 


0.001334 
0.001419 
0.001533 
0.001688 
0.001900 
0.002189 
0.002580 
0.003101 
0.003782 
0.004654 
6.005750 


0.575991 
0.575758 
0.575413 
0.574866 
0.574424 


'O 


h 


0.000988 
0.000970 
0.000972 
0, 000976 
0.000981 
0.000989 
0.000999 
0.001013 
0.001031 
0.001053 
0.001080 


0.118569 
0.118669 
0.118726 
0.118744 
0.118726 


We 4 


& 





0.000968 
0. 000390 
0.001014 
9.001041 
9.001072 
0.001109 
9.001158 
0.001207 
9.001273 
0,001354 
0.001450 


-0, 054616 
-0. 055247 
-0, 0553871 
-0. 056488 
-0. 057098 


- 
ow 


0. 300. 


82, 748 
82.817 
82. 870 
82.910 
82.936 





20 


Table V- 


Pb 
Eh 


0, 

0.500 
1, 000 
1, 500 
2. 000 
2. 900 
3.900 
3. 500 
4, 000 
4, 500 
5. 000 


29. 845 
29, 854 
23,360 
29, 262 
29, 861 





0. 887324 
0.959922 
1.038418 
1.122226 
1, 210522 
1, 302331 
1, 396643 
1, 492512° 
1, 589121 
1. 685807 
1. 782061 


4, 216877 
4, 206877 
4.196877 
4. 186877 


4.176877 


Normal pressure, p*24.50 Eh 


0, 017832 
0. 018778 
0, 019949 
0, 021376 
0. 023084 
0. 025090 
0. 027398 
0, 030004 
0, 032895 
0, 036054 
0. 039465 


—Q. 046 562 
—0, 048771 
~6, 048674 
~Q, 04727 
- 0, 046596 


4144 





0.017832 
0.019579 
0.021710 
0, 024302 
0, 027434 
0.031179 
0, 035603 
0. 040762 
0, 046704 
0. 053465 
0. 061077 


1.643312 
1, 647968 
1, 652300 
1, 656320 
1, 660081 


w 
$,.5 


h 





0.002001 
0, 002157 
0, 002359 
0. 002619 
0. 002956 
0, 003389 
0. 003942 
0. 004640 
0, 005510 
0, 006581 
0, 007881 


0. 587627 
0. 597366 
0. 596987 
0, 536502 
0. 595$17 


..” « 
' % 
bh 


0, 001316 
0, 001322 
0. 001326 
0, 001332 
0. 001340 
0, 001350 
0, 001364 
0, 001381 
0, 001401 
0, 001426 
0, 001456 


0, 1255388 
0, 125629 
0. 125677 
0, 425633 
0.125651 


0.001318 
0, 001349 
0. 001383 
0. 001421 
0.001465 
0, 001516 
0.001575 
0. 001645 
0.001727 
0, 001822 
0. 001934 


—0, 055255 
~0, 055903 
~0. 789545 
- 0, 057182 
-0,0570; 


Values of deflection coefficients for various values of axial compressive load in 
the x-direction, P, for simply supported rectangular plate, a/b=1.00, »=0.300. 


Cb 


0.975 
1,641 
2. 336 
3. 060 
3. 817 
4.604 
3. 421 
6, 268 
7.141 
8. 040 
8. 962 


85, $49 
85,923 
35,.$61 
86. 024 
86. 052 





al 


Table VI 


30. 436 
30. 446 
30. 453 
30.455 
30, 454 


- Values of deflection coefficients for various values of axiai compressive load in the 


x-direction, P, for simply supported rectangular plate, a/b = 1.00, :: = 0.300. 
Normal pressure, p = 30.00 Eh*/b*. 


Ww 





1.020802 
1.093572 
1.170845 
1.252071 
1. 336590 
1. 423682 
1.512635 
1.602790 
1.693571 
1.784501 
1.875199 


4.263910 
4, 253910 
4.243909 
4.233909 
4.223909 


“1,3 
h 


0.022490 
0.023659 
0.025052 
0.026690 
0.028587 
0.030751 
0.033183 
9.035877 
0.038824 
0.042011 
0.045424 


~0.060373 
~0.059704 
-0. 058926 
-0.058047 
~0.057077 


3.1 
h 


0.022490 
0, 024643 
0.027203 
0.030234 
0.033798 
0.037955 
0.042753 
0.048255 
0.054486 
0.061486 
0.069285 


1. 679468 
1.684135 
1.688463 
1.692479 
1. 696205 


“33 
h 





0.002652 
0,002878 
0.003161 
0.003516 
0.003961 
0.004514 
0.005199 
0, 006038 
0.007057 
0. 008282 
0.009742 


0.615745 
0.615448 
0.615031 
0.614503 
0.613874 


0.001515 
0.001620 
0.001626 
0.001633 
0.001643 
0.091656 
0.001671 
0.001690 
0.001713 
0.001740 
0.001772 


0.131466 
0. 131348 
0.131583 
6, 131574 
0.131526 


0.001615 
0.001654 
0.001696 
0.001744 
0.001798 
0.061859 
0.061930 
0. 002G61% 
0.062105 
0.002213 
0.002336 


-0, 055621 
-0. 056283 
-0. 056940 
-0.057594 
-0.058245 


eb 


1,291 
1,983 
2. 700 
3. 445 
4,217 
9.018 
2. 844 
6. 697 
7. 574 
8.473 
9. 395 


88.516 
88.593 
88.653 
88. 697 
88.725 





Table VII - Values of deflection coefficients for various values of axial 
compressive load in the x-direction, P, for simply supported 
rectangular plate, a/b = 4.00, » = 0.300. Normal pressure, 
p = 2.50 Eh’ /b*. 


Pb of We 4 Ws 4 a eb” 
Eh? h h h h 
0. 0.383 0.068 0.015 0.004 0.015 
0. 500 0. 395 0,077 0.018 0.005 0.516 
1.000 0. 407 0, 088 0.021 0.005 1.019 
1, 500 0.421 0,103 0. 026 0.007 1,522 
2.000 0. 436 0.123 0.034 0.008 2,027 
2. 500 0. 454 0.150 0.045 0.010 2.535 
3.000 0.473 0.185 0,064 0.014 3.049 
3. 500 0. 493 0. 230 0.099 0.019 3.575 
4,000 0. 508 0. 263 0. 188 0.020 4,137 
4.101 0. 498 0.256 0.238 0.015 4.276 
4.133 0. 488 0,249 0.264 0.011 4, 329 
4.151 0.478 0, 243 0. 285 0.008 4. 366 
4.163 0. 468 0, 237 0. 303 0.006 4, 396 
4.172 0. 458 0.232 0. 319 0.003 4. 422 
4.179 0. 448 0.227 0. 335 0,001 4. 446 
4,185 0, 438 0, 221 0. 349 -0.001 4, 469 
4.190 0. 428 0.216 0. 364 -0, 003 4. 492 
4.195 0. 418 0,212 0. 378 -0,005 4,514 
4.199 0. 408 0.207 0, 391 -0.006 4. 837 
4.204 0. 398 0.202 0. 405 -0. 008 4.561 
4.209 0. 388 0,197 0.419 -0. 009 4. 586 
4.215 0. 378 0. 192 0. 432 -0.011 4.612 
4.221 0. 368 0.188 0. 446 0.012 4. 640 
4, 228 0. 358 0.183 0. 460 -0.018 4. 670 
4.236 0. 348 0.178 0.474 -0.015 4.701 
4.245 0. 338 0.174 0,488 -0.016 4.736 
4.256 0. 328 0. 169 0. 503 -0.017 4.772 
4. 267 0. 318 0, 164 0.518 -0.018 4.812 
4.280 0, 308 0. 160 0. 533 -0.019 4. 854 
4.2094 0. 298 0.155 0. 549 -0.019 4.901 
4. 310 0. 288 0. 150 0. 568 -0.020 4.951 
4, 329 0.278 0. 146 0. 583 -0.021 5,006 
4. 349 0. 268 0.142 0. 601 -0. 021 5,065 
4. 372 0.258 0. 136 0.619 -0.022 5, 130 
4.397 0.248 0.131 0. 638 -0.023 5, 202 
4° 426 0. 238 0.127 0. 659 -0.023 5, 280 
4.458 0.228 0.122 0. 680 -0.023 5, 366 
4.494 0.218 0.117 0.703 -0.024 5. 461 
4. 634 0. 208 0.113 0.727 ~0,024 5. 567 


t 
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Table VIII - Values of deflection coefficients for various values of axial 
compressive load in the x-direction, P, for simply supported 
rectangular plate, a/b = 4.00, u = 0.300. Normal pressure, 
p = 7.50 Eh’ /p*. 

Pb wi, 4 “34 ef Ei 24 @b* 
Eh® h A Se ’ h? 

OQ. 0.961 0.223 0.064 0.918 0.115 

0.500 0.986 0.242 0.073 0.921 0.627 

1,000 1.014 0.264 0.085 0.24 1.143 

1, 500 1.044 0.288 0.099 G.929 1.663 

2.000 1.076 0.314 0.116 0.935 2.188 

2. 500 1,112 0.342 0.137 0.043 2.719 

3.000 1.15% 0.372 0.162 6.953 3.258 

3. 500 1.193 0.401 0.190 0.066 3. &07 

4,000 1.237 0. 430 0.223 0,981 4, 366 

4,500 1,282 0.457 0.259 0.100 4.938 

5. 000 i. 328 0, 482 0.298 6,122 5. 524 

6.727 1.416 0.499 0.469 0.250 7.716 

6. 780 1, 406 0.491 0,482 0.257 7. 796 

6. 805 1. 396 0.484 0.492 0.259 7. 838 

6,815 1. 386 0.477 0.501 0,259 7, 857 

6.815 1.376 0.471 0.509 0.256 7.863 








Table IX - 


1,500 


4, 000 
4, 500 
5. 00 


8. 6698 
8. 675 
8.677 
8.677 
8.675 


Values of deflection coefficients for various values of 


axial compressive load inthe x-direction, P, for simply 


supported rectangular plate, a/b = 4.00, u# 0. 300. 


Normal pressure, p #12.50 Eh Io*, 


Ww 
Li 


h 


1,373 
1,404 
1,437 
1, 471 
1,508 
1, 546 
1, 537 
1,629 
1,672 
1,717 
1. 761 


1,900 
1. 890 
1, 880 
1, 870 
1, 860 


. 352 
~ 344 
. 396 
- 420 
. 445 
- 470 
« 495 
. 921 
- 046 
oT 
- 595 


. 659 
655 
. 651 
648 
644 


5, 1 


118 
~132 
o 147 
. 163 
182 
202 
2225 
- 248 
274 
. 300 
327 


- 488 
- 486 
- 484 
. 483 
e 481 





. 037 
. 042 
, 049 
. 056 
, 985 
O75 
OT 
~ L00 
116 
134 
154 


- 491 
- 498 
. 594 
oli 
O17 


gid 


. 263 

. 189 
1, 318 
1, 8&2 
2, 392 
2.937 
3. 490 
4,059 
4.617 
5.194 
5. 179 


10, 617 
10, 640 
10,658 
10.672 
10,684 





Table X - Values of deflection coefficients for various values of axial 
compressive load in the x-direction, P, for simply supported 
rectangular plate, a/b = 4.00, u = 0.300. Normal pressure, 
p = 18.00 Eh’ /b°. 


Pb Wii "3.5 We 4 Wa. eb" 
or Se Sr 
02. 1,723 0. 464 0.171 0.059 0.448 
0. 500 1.755 0. 485 0.186 0. 066 0.984 
1. 000 1. 789 0.507 0. 202 0.074 1, 524 
1. 500 1.824 0. 529 0.219 0.082 2.068 
2.000 1. 861 0. 552 0. 237 0.092 2.618 
2. 500 1, 899 0.575 0. 256 0.103 3.173 
3. 000 1.938 0.597 0.276 0.115 3.733 
3. 500 1.978 0. 620 0. 296 0.129 4. 300 
4.000 2.018 0.643 0. 318 0.144 4.873 
4, 500 2.059 0. 665 0. 340 0. 160 5. 452 
5. 000 2.101 0. 687 0. 362 0.178 6.039 
10. 247 2. 305 0.795 0.544 0.578 12. 930 
10.254 2.295 0. 792 0. 543 0.585 12. 956 
10. 258 2.285 0.789 0.541 0.591 12.977 
10. 260 2.275 0.785 0. 539 0. 597 12.995 
10. 258 2. 265 0. 782 0. 536 0. 603 13. 009 


Table XI - Values of deflection coefficients for various values of axial 
compressive load in the x-direction, P, for simply supported 
rectangular plate, a/b=4.00, u = 0.300. Normal pressure, 
p = 24.50 Eh’ /b’. 

Pb 7) i "* i Ws Wo gb" 
Fh h 3 . . a" 
0. 2.053 0.569 0.224 0.083 0.672 
0. 500 2.085 0, 589 0.238 0.091 1.216 
i. 000 2.118 0.610 0.254 0.099 1.765 
i. 500 2.153 0. 630 0.270 0.109 2.318 
2.000 2.188 0.651 0.287 0.119 a. 875 
2. 560 2, 224 0.671 0. 304 0.130 6. #35 
3,000 2.261 0.692 0.322 0.142 4,0C3 
3. 500 2,298 0.713 0.341 0.155 4.574 
4,000 2.336 0.733 0. 360 0.169 5. 156 
4.500 2,374 0.753 0. 379 0.184 5. 732 
5.000 2.412 0.778 0. 398 0.199 6. 319 
11.964 2.675 0.916 0. 607 0. 659 15. 448 
11.973 2.665 0.913 0.605 0. 665 15. 478 
11.978 2.655 0.910 9. 603 0.672 15, 503 
11.981 2.645 0.907 0.601 6.678 15. 523 
11.981 2.635 0.904 0. 999 0.683 15. 540 











Table XII - Values of deflection coefficients for various values of axial 

compressive load in the x-direction, P, for simply supported 
rectangular plate, a/b = 4.00, u = 0.300. Normal pressure, 
p = 30.00 Eh’/b’. 

Pb “hd ¥3,1 “5, 1 “7,1 eb" 
Eh° . 7 . he 
0. 2. 286 0.643 0.261 0.101 0. 860 
0. 500 2.318 0.662 0.275 0.110 1.409 
1. 000 2. 350 0.681 0.290 0.118 1.963 
1. 500 2, 384 0.701 0. 306 0.128 2.520 
2.000 2.418 0.720 0.322 0.138 3.082 
2. 500 2.452 0.739 0. 338 0.149 , 3.647 
3.000 2.487 0. 759 0.355 0.161 4,217 
3. 500 2.522 0.778 0.372 0.1738 4.791 
4.000 2.558 0.797 0.390 0.186 5. 369 
4, 500 2. 594> 0.816 0.407 0.200 5.952 
5.000 2.630 0.835 0.425 0,215 6. 539 
13. 331 2.920 0.996 0.651 0.724 17. 476 
13. 337 2.910 0,993 0,649 0.730 17. 502 
13. 339 2,900 0.990 0.648 0.736 17, 524 
13. 339 2,890 0.987 0.646 0.742 17, 542 
13, 337 2, 880 0.984 0.643 0.747 17. 556 
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Table XIII - 


Values of (Pb/ Eh°) critical at various values of a/b 
and pb*/Eh*, These are results using equation (8). 
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Table XIV - Values of (Pb/Eh°) __ 
Critical 


and pb /Eh*. These are results using equation (7). 


ae Various valmesvor a/b 





eine eon = | 

4. 8. 68 10,26 11.98 13.34 
3. 8.77 10. 33 12. 03 13. 37 
3, 8.87 10.41 12,08 13. 40 
3, 8.97 10, 50 12.14 13.45 
3, 9.08 10, 58 12,21 13. 50 
3. 4. 9.64 11.04 12.57 13.77 
3. 4. 10.00 11. 32 12. 80 13.96 
3. 3, 10. 00 11.65 13.06 14.19 
3. 3. ae 12. 50 13.74 14.77 
3. 3. 8.98 13.09 14.17 15.13 
3. 3. 8.97 12.47 14.26 15.21 
3 3. 8.96 11. 69 14.47 15. 38 
3. 3. 8.96 11. 42 14.70 15.57 
2. 3. 8.98 11.19 15.03 15. 82 
2. 3, 9.03 11.05 13.98 16.11 
2. 3. 9.10 10.97 13. 40 16. 47 
2. 3. 9.19 10. 93 13.08 -| 15.19 
2, 4. 9.79 eee: 12.71 13.98 
2. 4. 9.91 11.25 12.74 13.95 
2. 4. 10. 24 ikea 12. 87 13.97 
2. 4. 10, 83 11.97 13. 20 14.18 
2. 4. 11.58 12. 60 13.72 14.61 
2. 4. S415 12.99 14.06 14.90 
2. 4. 8.81 13.44 14.45 15.26 
2. 4. 8. 64 11.64 14.91 15.67 
2. 4. 8.56 10.93 15. 43 16.16 
2 4. 8. 56 10.61 13.12 ee 
2. 4. 8.61 10.45 12.55 14, 30 

8. 69 10.45 12.15 Sees 


10. 06 11. 36 12.71 ih) 13.74 
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V. DISCUSSION OF RESULTS 


The most remarkable part of the results was the way tne 
critical buckling load behaved with respect to varying values of 
a/b (Fig. XXXIV). A= can be seen in the corrolary figure (page 64) 
to Fig. XXXIV, there is a striking resemblance to B. yan's 
buckling load solution of rectangular plates solely under the 
action of edge loadings. It may be recalled that the wave-like 
pattern of B: yan's curve of K versus a/b was due to the nature of 
the K-formula anidto the minimizing procedure adopted in obtaining 
the critical load for any given value of a/b. Bryan's K-formula is 
a function of the plate's dimensions and of the number of half-sine 
waves the plete takes after buckling. Thus, it would seem from the 
foregoing statements to presuppose that the patterns obtained in this 
thesis are also due in part to the number of half-sine waves. 


At first glance the results would give no indication as to the 
number of waves the plate takes at the incipience of buckling. But, 
checks had been made by the authors on the relationships between 
deflections at the midwidth of the plate and at various distances along 
its length. In other words plots were made of the deflection ratio, 
w/h, versus the distance ratio, x/a, at y = b/2. These were done 
for different (a/b)'s at the same normal pressure. For example, it 
was found that for a p = 18.00 Fh* fte* there were five buckles at 
a/b = 2.40. (The two half-sine waves resulting from the initial 
general downward deflection of the plate due to normal pressure 
were not inclu.jed in the count.) When these checks were expande to 
a wide range of values of a/b, it was found that points lying on the 


same trough of a p-curve gave the same number of buckles. 


Since only representative tables are given (Tables I to XII) it 
was deem-d unnecessary to insert all these checks inthis thesis. If 
verifications are to be made on this particular point the authors suggest 


going through the programs shown in Appendix B. 
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There is, howeve:, no positive explanation as to why only 
an odd number of buckles is possible. There are either one, three, 
five, etc., buckles. An even number of buckles, which are obtainable 
in the case of Bryan's, seems entirely out of the question. On this 
particular aspect, a possible explanation that the authors believed 
could have happened was in the manner the deflection equations were 
obtained. The approximated deflection equations were follow-ups of 
those use in the references | 2}, [3], [4]. The equations only have 
the odd-numbered, subscripted deflection coefficients which, in the 
final analysis, resulted in an odd function for the w-equation. At this 
point, there is no way of predicting what the results will be if the even- 
numbered subscripts of the deflection coefficients are also taken into 
consideration. 


In reference to the reliability of the results, it is to be noted 
that the basic theory is vali only within the range of Hooke's law. 
No allowance was made in the original differential equations on the 
change of Young's modulus when the axial load exceeds the propor- 
tional limit. Hence, all the results of this thesis presupposed that 
the material had been within the elastic region at all times. Experi- 
mental results are not available on the case of plates with the same 
loading and boundary conditions such that, definite identifications 
cannot be made as to whether the results are within the elastic 
region or not. 


When a further study is made of Fig. XXXIV for a/b less than 
Or equal to one, there is a very wide discrepancy between the calcu- 
lated critical loads and Bryan's that it is quite doubtful whether the 
expected critical loads are below the yield strength of steel for a 
majority of plate characteristics. 


The following computed results are hoped to illustrate the point: 


Consider a plate with a b/h = 104.7 and acted on by a pressure, 
p= 12.50 Fh*/p4 = 3.11 psi. The b/h may correspond to a 14. 02- 
pound plate with b = 36 inches, andthe p may correspond to a hy4dro- 
static head of 7 feet. 
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Case 1 Case 2 Case 3 
a/b = 1.50 a/b = 1.00 a/b = 0.60 
K = 12.8 K # 28 K = 78 
C. -- $1800 psi o. = 69000 psi Go. = 192000 psi 


If mild steel or high-tensile steel (HTS) were used with yield 
stresses of 33000 and 47000 psi, respectively, it is apparent that 
cases 2 and 3 exceed their corresponding yield strengths. Hence, 
the theory feils for these cases. If HY80 or HY100 were used then 
theory is found to fail only fur case 3. At any rate, the fact remains 
that the basic theory does not differentiate between the elastic and 
inelastic regions. A consequent conclusion, therefore, is the need 
for a check on the yield strength whenever the formulated o. of this 
thegi; is used. 


Except for the restrictions mentioned in the previcus paragraph, 
it is obvious that the results obtained confirmed the conclusion of Levy 
on the effect of normal pressure on the buckling load. But, it is also 
evident that to neglect the normal pressure in predicting buckling strength 
is unwise design from an economic point of view. 


Lastly, the computer programs are limited in scope, although 
they worked properly in providing all the needed results of this thesis. 
They were the consequent developments of the approximated deflections 
equations. A more desirable program would be to start with the 
equations shown in reference [ 2]. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 


1. Normal pressure always increases the buckling strength 
of rectangular flat plates. However, the buckling loadsa obtained by 
means of Levy's golution are not conclusive values; that is, the 
buckling load may or may not be within the proportional limit. In 
other words, theory does not differentiate between the elastic and 
inelastic regions. 


2. Plotted values of critical loads against the ratio, a/b, 
for a constant value of normal pressure exhibits the same trend 
and characteristics as the plot of Bryan's formula for zero 
pressure. 


3. To neglect the effect of normal pressure in the determination 
of buckling load is unwise design from an economic standpoint. 


4, In general, the buckling load formula can be written in the 
form, 
2 2 


aoe = 


© 1844-0") 


where, this time, K is not only # function of plete geometry but also 


o 


a function of the normal pressure and Poisson's ratio as well. 


Fecommendations 


1. A better representative plot of K versus a/b is needed by 
including the even-numbered subscripted deflection coefficients of 
Levy's solution. 


2. In consonance with the first recommendation, a :evigsed 
computer prcgram is Jesirable, possibly to start right at the series 
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forms of the equations shown in reference {2}. This step would 
eliminate using approximate deflection equations. 


3. It is recommended that a study be made of the most 
feasible experimental set-up to simulate the boundary conditions 
stated in this thesis; and hence, the ensuing experimental 


verifications. 
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APPENDIX A 
DEFLECTION FQUATIONS 


The six (6) simultaneous cubic equations involved in the solution 
of the deflection equation, 


Ma -, WX. Ty xX omy 
Ww Wid sin > sin 5 77 5 sin i sin b 
37x HEX — 
+ We 4 sin ‘- sin +++ We 3 sin ———sin 


mx ,., Sy omx 1 
+ W, 5 Sin 4 sin 5 + we , sin ; sin 


are as follows”: 





a eek ow 4 es 16 Pb | 
6 Eh at Nt 6 Fh “1,1 
T B(i-u") “3 

3 we 2 

+ (a, + ao) Why - 8a, w 1,1 %1.3 7 38, Wy 3.1 

+ Qa, + 4a. + 16a, + a 3) Wy 1 2 3 + (4a,+ 9a, . 16a. 
2 2 

+ a4) wy 13. i* 32a, Ww 1,113 W3 4 de 81(a, +a.) Wy W3. 3 


- 18(a, > ag) w - 18(a, +a,) Ww 


1,1%1,3 %3.3 1,1 %3,1 “3, 3 


‘ 2 
~ (9a, + 64a, + 25a,)w) owe , - (9a,+ 64a, + 250,) “1,676.4 


+ (36a. + 36a, - 225a, +225a,) Ww 


, ~ 6(a, 


1,3%3,1 “8.8 


* 3a9) W113 Wy 5 ~ Gla, +3a,) “4,1 "he 76,4 


+ 162a + 162a, Ww + (25a, +49 


61,1 “s,3 1,6 1,1 “S878, 1 2 


*The corresponding significance of each suoscripted a is shown at the end of 
this Appendix. 
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256a, + 49a + (Za, + 64a, + 2562 
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2 ae Pk 
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PoOeeTTON OF THE DEFLECTION EQUATION 
We Wisin sin + W5,4 57 = sq as 
+ Ws, sin 2X ca sim Je “2 + Wz) sm stall sin - 
hacks <a fo SE Orr rer RA CONTINUOUS 
GAeNGe IR @weATALC COAD =ROM ZERO TO A VALUE @ITHIN 
ifesvroernar OF THE CRITICAL BUCKLING LOAW. ZERO 


PemeRAL OAD 1S! EXCLUDEN. 
MAIN PROGRAM 


XEQ 

LABEL 

larS 16 
FiQiReeRuan 

READ 1s IND 


i FORMAT (15) 


Pron T 2 
2 FORMAT (82H1L PNORM AXIAL WCls1l)/H W(391)/H "Ss 7 ee 
m7 8 6STRA EN ERROR ) 


DO 44 JOB = 1s5:ND 
READ S3°U, ABPHAS PNORM, LOOPs LAP 
3 FORMAT (3FC8 88 212) 


DIMENSION 81(4)+ B2(4)5 B3(4) 





CALL CONST (Use ALPHAs PNORMs Als A2s A3s A4s Ads ABs Als ABs ADs 
1 A1Os All>s \12s A1l3s Al4s AWSs AlGemeey > Bie Seer 


C PNITIAL ESTIMATES OF Aa 
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DIMENSION W(4) 
W(1) = 04400 


WC2) 0.300 


W(3) = 0.200 


W(4) 


02190 
THE SOLLOWING IS THE ROUTINE FOR 

GEUCUCMWTNG THE BEREECTTION GweFFICIENTS FOR 
CHOSEN VALUES OF AXIAL LOAD AND NORMAL LOAD 

DIMENSTON A(4.4)5 Cl4)s Xl4)o 26305 $1354) 

AX = O« 

MG = 1 

96 WO OLL = 153 

CALL SESAME 

GME CROW TO IIEC sx dM) 

GO TO (55 42)5 M 


ERROR 


Oe 


DO 6 J = 194 


Breen = ERR@R+ ABSF(X(J)) 
Meh tERROR -M7000]) 9: 95 7 
DO @ T = 154 


wel) = WCl xT) 

GO TO 4 

CALL CPRNT (PNORMs AXs Ws ERRORs Al) 
Gee SVORE Wes Ws Zs Ss LL) 

AX = AX + GO85 

CONTINUE 

DEL = Oe5 


KO = 1 


: 
> . 
ae 
_ 
a e 
=A ‘ 
af 
» Cin —_ 





a? “ae TO (239 255 27)» KO 
eo Wels 0805 

meres 1. 2'154 

we ay = S(2391) 
24 $(33s1) = WT) 

px = 702) 

Z(3) = AX 

AX = AX + DEL 


KO 


2 
KK = 9 
MmeLlL VeXCH ( Z9er) 
ce TO 13 

m> DEL = 0.01 


DO 26 IT = 194 


wel) a= S291) 
26 S$(391) = W(T) 

AX = Z(2) 

Z(3) = AX 


AX = AX + DEL 


KO 3) 
KK = O 
Cm, CMH (2s 3S) 
Go *O 13 
we “POLLOWING IS THE ROUTING et Se 


PeMGINITY OF THE BUCKLING LOAD. Wl) Si Sere 





THE INDEPENDENT VARIABLE IN PLACE OF SIRE Axia 
. LOAD. 


ma DO 28 “1T = T54 





ERROR = 49. 

DO 32 J = 194 

FRROR = ERROR + ARSF(X(J)) 

IF (ERROR - 04-0001) 36:2 36:5 33 

me = KK + id 

IF (KK - LOOP) 34 5 34s 40 

pemes it. = 256 

wary = WT) + ACI) 

pee = AX + X(1) 

GO» T O30 

CALL CPRNT (PNOPMs AXs We ERRORs Al) 
IF (AX - 2(2)) “%4s 445 37 

mes) = AX 

DO 38 I = 194 

S$(391T) = W(T) 

Ore ie 7(5) -— AZ) 

CONTINUE 

GO TO 44 

PRINT 41 

FORMAT (40H INFINITE LOOP OR NEEDS MORE ITERATION) 
CALL CPRNT (PNORMs AXs Ws ERROR» Al) 
GO TO 44 

PRINT 43 

FORMAT (19H @rOouT WENT CRAZY) 


CONTINUE 





Gyre EXIT 
COMMON Ws AXe Bl» B2s B3s Als A2s A@ts AGs ASewetoe ARe ACe AlDs 


] Ads Al?s AlAs Al4s AZSs ABS) s\) eee C 
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END 
SUBROUTINES 
Supreme ot ine FOR SOLVING THE LINFARIZED SYSTSsM 
OP PEQUATIONS BY THE CROUT REDOX METHOD 
BABEL 
EgSTS 
FORTRAN 


SeSROUMINE CROUT (As Cs Xs M) 

DIMENSION AlGst)s Cl4)s AA(494)5 CC(4)s X(4) 
DO 1C0 I = 194 

AA(Ts1) = AlTs1) 

WO 191 Js 254 

AA(1lsJ) = AllsJ)/4A( 151) 

DO 1°02 TI = 294 

DOere2 J = 258 


fea ese = ON 


B® 405 I = 1134 
lt = J - 1 


DO 104 K = 1sLIMl 


Beets J) = Fehrs J) + (AAC IK) * AACK9J)) 
AAs) = Mie) — AA(I.J) 
Le Pra ( J 1065 1163s 106 


ae ~ Joo 9 10s 107 


a) 


YJ ‘ ' 


75 = | on fet 
ae! 








—- + 1 

Jor=* J 

POTVTO9o tJ = JIs4 

LIM2 = sia 

Me wos K = 1sliv2 

AA(TsJ) = AAC(IsJ) + (AACTsK) * AAC K9J)) 
Met l>J) = Alisa Aare 9 
Rat's J) = AAUT SO) AAT) 

Jel + 1 

mo TO 103 

we le 6€6hdT 6 = (U1 4 

ml) MPO. 

Sery 7s wd, 

eer 2ec(1)78AC1 > 1) 

DO 113 I = 294 

LIM3 = - 1 

POP ll? eK = bok MS 

moe CECT weeAA( 1s) ™® CC(K)) 
Sewer] CC eee cr) 

Gort) = CCR) 7AATI >I) 

X(4) = CC(4) 

WeeilS It =F1s3 

II = 4 I 

LIM4 = II + 1 

DO 114 K = LIM4s4 

ee = Xie (AAC TI sK) * ATK) 
oe tye X(T!) 


ic ia | 
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GGnmO. 118 
PRINT 117 


FORMAT (28H SINGUEARTTY = NO SOLUTION) 


RETURN 

END 
SUBROUTINE FOR CALCULATING THE NON-DIMENSICNAL 
STRAIN AND THEN PRINTING ALL THE NEECED RESULTS 

LABEL 

Bas T 8 

FORTRAN 

SUBROUTINE CPRNT (PNORMs AXs Ws ERROPR®s Al) 

DIMFNSTON W(4) 

SUM = Oc 

BO 70@" IJ = 134 

Meee ces IJ = | 

SUM'="'SUM + (XY * W( TJ) )*#2 

STRAIN = AX 4+ 341416*%*2 * SUM/(8. * Al) 

PRINT 2012 PNORMs AXs (W(T)> [=194)9 STRAIN>s ERROR 


FORMAT (F8e2s Flle5s 4F10035 F1lle5s F125) 





RETURN 
END 
Smee OUTTNE FOR STORAGE @F TAREE StlseOeyvVebue. 
OF AXIAL LOAD AND THF FOUR DEFEEC TIONS CCEr ler. 
LABEL 
Basile 
FORTRAN 


SemROUlTINGe=STORF (AXs Ws Zs Ss LL} 








DIMENSION W(4)5 Z2(3)5 S(3594) 

Zep R= AX 

00 203 I =154 

SL Lel }e=oWCT) 

RETURN 

END 
THIS SUBROUTINE IS JUST A SHIFTING PROCEDURE 
GF VALUES FROM ONE STORAGE LOCATION TO AWOTHER, 

LABEL 

e1sne 

FORTRAN 

PUOBRGUTINE EXCH {(Z; S) 

Seve mS TON 263)6 S( 394) 

PO ame €6T = 154 


Der 704 J =~ 1592 


wos!) eg 3( JSideT) 
Die O25 41 “ele 2 
J@e= IJ + 1 
wer.) = 7 0 
Deas JST Seis 
Saco JI) = 
Z2G3) = O- 
RETURN 
END 
SUBROUTINE FOR COMPUTING ALL THE NECESSARY 


CONSTANTS IN THE MAIN PPOGRAM, 





| PASE 





res 





FORTRAN 

SUBROUTINE CONST (Us ALPHAs PNORMs Als A2s A3s 4s Ads 459 ATs ABs 
1A9s A1Qs Adds A122 Al3% Al4> Al5® \169 Al7>s Bl» 829 92) 
DIMENSION B1(4)»s B2(4)5 B3(4)s P(9) 

P(1) = ALPHA#*2 

exe y = 1le/P(1) 

Der 3ago ft = Tbs7 

Yl = [*#*2 

Pl <= Seeee2 

BOO PWT) = P(1)/0P01) + Y1)**2 

or FPP 1) 

A2 = P(l) + P(2) 

ose) ) Hee. PD) 

a= Piel) + 695. ¥ P(2) 

A> = P{1) + 24lew * P(2) 

oe FP 64, PE 3) 

Roe ees * Pe 76. * (3) 

ae = 30's PTT) + Gree ea | 

CS Sarre + 256 HOG) 

AYO £935 *# PP) + 64. * P(5) 

M11 fobs FeRH]) e+ ee ey + 146 * (3) + P(4) 
ieee eee + 2. YP (2) + 2EHe * PE) + :eRGD 

AiR = Ge +e) + 1epee * PI?) + 1206e * PC 2) FPR) 

A164 = Ge 5a) + 250 * P(2) + Ble * P(4) + 160 * P(5) 
ATS Ge * P(1) AGL. * PU?) 6254 #9 Pe 16.6 P(7) 
Al6 4e * P(1) 49, * P(2) 2560 * P(5) Bl. P(6) 
A17 bo * PCI) 14¢, * PC 225, * P(4) oO. P(6) 


@ _ 
'F 
) a e 
— aE” == 
‘ 
$ = 
= 
. 5 
; | 





—s 
ty . 
a= 
ie ‘One 
wr et <8) pe ee = Tee 


me A) el 



























io $n) Pee ene 


'» , 1 4 














‘i 








MO “3a 1 I = 154 
2 2=@2* |] we] 
Y2a9e 2 * [°# 91 


ml (.) 2564 * P(1) * PNORM/S(Y?2 * 32.41416%%6) 


a2 (1) Ame / OS Cyr erUuseD) FERC 127) 


BO? B3(1) = =(3OG'# Y2#%2)/(821416**2 ) 

RETURN 

FND 
SUSr OUT] NEeFORS COMPUTING THE COEFFICTENTS 
One tino SS OF Thm CROUT REDOX METHOD 

Lm EL 

BISiés 

FORTRAN 


SUBROUTINE SESAME 
COMMON Ws AXs Bls B2s B32 Als A2s A3s A4&s A5s Afy A? ABs AGe AlDs 
} Alls Al2s Al3s A1l4%s A155 Al6s A17s MGs As € 
DIMENSION W(4)3°B1(4)5'°B2(4)5 B3(4)5 Al4s4)5 C4) 
GO TO (401s 400): MG 
400 A(191)=B3(1)#W(1) 
Al29s1)=B3(2)#W(2) 
A(391)=83(3)*W(3) 
A(431)=B3(4)#W (4) 


GO TO 402 





403 Al191)=2=B2(1)4B39(7) *¥AX43, FAPEN(1 ERP —KA FAT EW( 1) NID) 4A1 1 FN (9) F¥D 
1 FAL GEW( 32) HHD4+ALBEHW( G4) FHD-F7 HAR EW(DIEW(2)-2 EF ALORV (2) HN( 4G) 
AC PST) =-3 oe ew (1) et A+2 , FA 1 PW SO Ane OT) es) 
i +2 o ¥ANBEW( 7) EN (2AS—2eHADEW (2) EN (4) +A] 7T*W (3) RWS) 


N( 251) a—PoeFAQHW(T)HW( DDE ZOHATGHW( TIEN (2)—7.eALTORW (1) HN (4) 


‘ ;e. ‘ 
e 


7 e *." 


7. =e . fs 


Sl - 














A 


1a 


Lim 
| ipeeiw owe ; 
- a 





ye | bel tel =m trees: we whtP® 


peered Loree s 
‘eae eee 


, 6A " ot='' 


7 a 
<= 


Ae f 
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1 +AG#¥W(2)¥#24A1 THEW(2)*W(G) 
A491 )=-2 o #ALO*RW(1) FWA) +7 HAL EAW(] )HW(4)-AQHW( 2) H¥D4A17#W (2) HW) 

902 Al192)=—-2e#ALEW(1) HH 242. 7ALLEW( 1) HWI9) +2 .*ABHW( 2) HW(3) 

1 —2e*AID*EW(2)FWI4)—2e*AB*W( 1) FW(3)+ALTEN (2) HW(G) 
A(2927)=B2(27)4B63(2)*¥AX+A1 LEW] )H4D4+2 XAQZ*W( 9D) ¥HD4+A1 DPEW( 3) F¥D 

1 +A15*#W(4)##D+7 oe #AGHW( 1) ¥WI3)-DeRADKW (1) WIG) 
N( 32927) =-AREW( 1 )#*94+2 .#AK RWI 1) EWI 2) 47 HAL PEW( 7) HW 2) 

1 +AL7TEWCLY#W(4)42,*AT7THW( 3) #W(G) 
A( 497 )=—P ae XFAQDHW( 1 )PHWID) 42 Ce KAITSHW( 7) AWG) +ATEW( 2) HXD4+4h1 THEW (1 YEW) 
A( 193 )=A6*#W(27) #247 -HALGEW( 1) #W0 22-2, #ARHW( 1) FWD) 

] =—2 e*¥ALO#XW( 1) #WI4)t+AL7TEW(2) RWG) 
A(2932)=-AS#EW (1) HHD47 HAI ZEW(2) W092) ADC KATH»>WN(3)*W(G) 

1 +72 o#AGKHW( 1) HWI2)4+A1 7THEW( 1) HWOG) 
A(3293)=BR2(2)+B2(2)#AX+ATGEW( 1) HH P+ALAHW( 2D) RE D+3CHAGHEW( 3) HRD 

1 +AL3Z#W(4)*#24+2 eX¥ATEW(2)¥W(G) 
A{493)=-A1O*¥W(] )F*FD47 FA TEW( 2) W020 4+90%ATZHW( 2) HW( G4) +ALL THEW) RAD) 
Al1s:4)=-A9#W (OD) #9470 .% ALT 6EW( 1) FW(I4)-9- FALORHW(] ) W039) 4A17#HW( 2) HH 3) 
A(294)=N7#W( 2) HHD47 CHATS EW( 7) WIA )-2. KADEN (1) FWD) 4+A17#W( 1) FW) 
A(394)=-A1O#W( 7) #42742 eFALAHW(2)HWKIG)FALTEW( 1) #W(2) +2 0 *ATHW( D2) HW( 2) 
A(494)=B2(4)4+B3(4) ¥AX+A] 6¥W( 1) H¥D+ATRHW(D) HHD+ALAHW(2) 4H] 

i +3 .%#A5SH#W(4) X*2 


C(1)2=R8161)—-B2 (7 HW01) -R9 01 HW) AXA DEW( 1) HH ALA HALTHEW( 1) HHDEW (9D) 





1 —AL1T*EW( 1) *WO7) HH D—AGHEW( 7) HK DREW (2) +AQKEW(D)KHDEW(G) 
? —A 1 GEWP eW 6 2) HE D—-A LA EW YEW (4) ERO +2 et ARM DN Oe) a) 
z) +P eoX*#ALO*W (1 )H#W( 3) ¥WI4)MKALTEW( 2) ¥W(A) RWG) 


C(27)=B1(67)—BR9(7)HWI72)—-RALD)PHAXHW(D)Y+FALHEW( 1) HHA2-ALTEW( 1) HDHD) 
1] tAREW( ] )HHDHW (2) -AQZ*W( 7) KR HA—ATDHW (9) EW (3 PRE D—ATHW( 2 ) eH DEG 


? —“AYGEW( DH CL) HEI—F HARRI (1 PHW(D )HWI 3) +7eFACHW( 1 HID) HAG) 










| / F thi tie ae 
et ae 
ma taneses(ie ee SEOUL a pe 


si were Sipaear BF oeriae 


it Ce " » af OT =) 
iow —_—i om ssealoewse! ' 


‘baer 


=> 












3 —AL7T#W(1)*W(3) #W04) 


C(2)=B1(3)—-B2(3)#W(32)—-B3 (3) ¥AX#W(2)4+ARHW( 1) HDHD) 


1 —~A1L4G#EW(1) #¥2*#W03)+ALOEWO] HED EW (4) -AE*W( 1) AWD) HD 
? —A1LP*W( 2D) HK APHW (32) KAGHEW(A) KHAKATAHW (2) HW( 4) KHD 
3 —ALTEW(1)#WO7)EW04) HD eHATHW( 2) HW 3) HW) 


C(4)=R1(4)-B2(4)4#W(4)-B2(4) #AXHW(4)+A1OFW SY HH DHW (3) 


1 —~A1L6*#W(1) F42RW(4)+A9RWO] AWD) HH 2-ALSHW( 2) HXDEW (4) 
2 —ATHEW(2)HW03) XH D—A1 ZBEW(3) FRDEW (4 )-ASHW(4) HR 

3 —AL7TEW(1 )¥W02) FW) 

RETURN 


END 








M3888-3401 sFMS »DFBUGs1351250020 


SOLUTION OF THE DEFLECTION EQUATION 
: UX sin U3 AX oy SL 
W = Wi.on— sm +Wy 35M —S sen =e 


31x é 

+Wzq sn 2E mAh. + V5.3 sin SE sim 210 

3% gq TH 
ab 


THIS METHOD 1S SFT™UP FOR A CONTINUOUS 


Tx SHY 
+ W5 sim 3 si) — +Ws | sin 
CHANGE IN AXTAL LOAN FROM ZERO TO A VALUE WITHIN 

THE VICTIM OPMTHE CRITICAL "BUCKUTNG LOWD. ZERO 


aero OAD 1S EXCLUDED. 
MAIN PROGRAM 


XEQ 
EAE L 
2ePs7s 
FORTRAN 


READ 1+ IND 


mOrerAT (15) 
Pra | 2 

2 FORMAT (114H1 PMNORM AXTAL W(121)/7H Wiles Seva WO351)7 
1H W(393)/H W(1s5)/H WO shee STRAIN FRROR ) 


DO 44 JOB = 1sIND 

READ 25 Us ALPHAs PNORM, LOOPs LAP 

PORMAT (3'F Garerc 18 ) 

DIMENS TON WOM, R454)5 Bl@)s B2(5)>+ B3(3)>2 AlGs6)5 CIiG)s Atay 
CALL CONST (Us ALPHAs PNORMs Ro Bls B2s B3) 


INITIAL ESTIMATES OF W(T)'S 


91 





W(1) = 902100 

W(2) = 02010 

Wt3) = 02910 

W(4) = 02901 

W(5) = 02001 

W(6) = 042001 

MeeerOlLOWING IS THE ROUTINE FOR 

moeoiereernG THE DEFLECTION COEFFICIENTS FOR 
CHOSEN VALUES OF AXIAL LOAD AND NORMAL LOAD 

PPMENSTOM 2(3)> S(356) 

AX = O-. 

mer = FL 

RONE = R(1) 

Corre SLL = 153 

GALL? SESAME 

eee GROUT (Ase Cs Xs M) 

SOTO U5 5s 42) 9 

ERROR = O« 


poo SC CU 


error = ERR@R + ABSF(X(J)) 

Br (ERR@R = OFD00001) 92 95 7 
Po™ | = 16 

mel yo= WUT) ee x(t) 

GO TO 4 


CALL CPRNT (PNORM, AXs Ws ERROR® PONF ) 





CAUCL STORE (AX 9 lf» Ls Sis Ele) 


AX = AX + 025 


BO CONTINUE 





KO 


iT 
ay 


KK = O 
B® 12 7 = 136 


fe] (eee ol(isl)y7(Z(2) - 2(1)) 


Bem Coot) = S(297))/712(3) - 2(2)) 


meee (202) = 21g) / 26 


F4 


(eae — F1)/Fe 

ao -ereery tere = DE 055.* F4 * DEL **2 
mali xe tzs S) 

CALL SESAME 

BREE GROUT (As Co Xo M} 

GO TO (145 42)>5 M 

ERROR = OO. 

i To ae = 156 

MenOe = ERROR + ABSF(X(J)) 

me .errROR — 0.000001) Peer Tos 6 

KK = KK + 1 

peer = LOGE, 17s 17s 40 

DO 18 %- = 196 

MerTy = WIT) + X(T) 

GO FO 13 

Me "OWL 1) -— S251)) 223 225 20 

CALL CPRNT (PNORM, AXs We ERROR: RONE) 


Z(3) = AX 





DO 21 I = 196 
me S(3e1) = W(T) 


AX = AX + DEL 





cha 0 

Ge Te il 

Ge yo (285258 27)+ KO 
DEL = 0605 

DO 24 I = 196 

WETY = S€251) 


S(3s1) = W(T) 


AX = Z(2) 
Z(3) = AX 

eX e= AX Fe°DEU 
KO = 2 

KK = 9 


Game NCH ( 2 sey 
Go TOP 13 

Pel = 0.201 

Cewze 1 = 156 
W(T) = S(2s1) 
S(3e1) = W(T) 


AX = Z2(2) 





Z(3) = AX 
Pex oe AX + DEL 


KO 2 


—_— =. 

CALL EXCH (Zs S) 

Go TO 13 
C THE FOLLOWING IS THE ROUTINE WHEN IN THE 
C VICINITY OF THE BUCKLING LOAD. W(1) IS MADE 


C THE INDEPENDENT VARIABLE IN PLACE OF THE AXIAL 





LOAD. 


27 ee 280 6d Ce OO 
Sue] ) = Sl2sI) 


S(2s1) = SlosTI) 


mes WUT) = S351T) 

mo! @ ete) 

mx2) = Z2(1) 

Am =) 2(3) 
ME FOLLOWING VALUES OF DL1 AND DL2 ARE ARSITRARY. 
Wreee ROT VALUE OF DL1 IS £O CHOSEN SUCH THAT @3een 
Pome emer nsel BY Di2ds THr CONSEQUENT I TEReiae 
SOLUTION WILL NOT GO BACK ON THE PREVIOUSLY COMPUTED 
PATHe eNCwre = THe WUT) OR THE FIRST HARMONIC STARGS 
WOSDECREASE WIN VALUE SEFORE THE CRITICAL BUCKLING 
LOAM IS REACHED) 

Pe IS= 04620 


eZ = Owl © 


~ MG 


N 
N) 


Wem> LAST DO LOOP IS HOPED 70 GO BEYOND THE CGRT- 
TICAL RANGEe START WITH LAP = 306 
WO'39 LAST = 1+LAP 


KK A 


Ao@e= AX + OU] 
W(1) = W(1) =- DL? 


DO 29 J = ZO 





29 W(T) WOT) + (S639) -— S29J)) 


CAREW cH (23s S) 


eo CALL SESAME 








4] 


4? 


43 


4k 


CARE CROUTECASECs Xs M) 

GO TO (315s 42)5 M 

ERROR = 0. 

wo 32 Y= 156 

Ele ORe= ERROR + ABSF(X(J)) 

[Pe CERReK - OS Gme001y 3657365 33 

he= KK + 1 

IF (KK - LOOP) 34s 345 40 

DO 35 I = 296 

W(T) = W(IT) + X(T) 

ax = AX + X( I 

GO» TGm3.0 

CALL CPRNT (PNORMs AXs We ERROR» RONE) 
IF (AX = 2(2)) 443 445 37 

Z(3) = AX 

DO 28 IT = 196 

Bpeowne = W(TR 

Pee - £'8) = Rene ) 

CONTINUE 

GO TO 44 

PRINT 41 

FORMAT (40H INFINITE LOOP OR NEEDS MORE ITERATION) 
CALL CPRNT (PNORMs AXs Ws ERRORs RONE) 
GO TO 44 

PRINT 43 

FORMAT (19H CROUT WENT CRAZY) 
CONTINUE 


CAGE eT T 





COMMON Ws AXs Bls B2s B39 Re MGs As C¢ 


END 
SUBROUTINES 
SUBROUTINE FOR SOLVING THE LINEARIZED SYSTEM 
OF PEQUAT Pen eS BY “THEPCROUT REDOX METHOD 
LABEL 
' EFS 8 
FORTRAN 


SWBROUTPNE CROUT (As Co Xo M) 
DIMENSION Al696)s C6)» AA(696)9 CC(6)s X(6) 
1m 86 = lb 


HOO WAI 51) A(T91) 


DOV TOT FY = 2456 


mOl MAC) sJ) 


Al@em@) 7A 151) 


Vere | = 256 





OOmeoe J = 256 


HO2 AA(IsJ) = O- 
J=2 
mOS Tit= J 


veer > [tert sé 


LIMI] = J 


' 
+—s 


DOO a wkK = T4LIM) 


meer AATTsJ) = AA@IsJ) + (AACIsK) * AA(KsJ)) 
me AA(Tsu) = AlTs@) — AAT >J) 
' eee A( J 5 Jom) 1065 116. 106 


mero IF (6 — J} ~ 9105 1105 107 








nO? fF ee J 
a= J +e) 
Jae J 
pO m09 Jae Jiwd 
LIM2 = =e 
Oe 106 K = 15LIM2 
mom AATI1sJ) = AA(TsJ) + (AAT sK) * AA(Ks9J)) 
Reet ss) & ACT #9 ABRs) 
109 AA(TsJ) = AAC TsJ)SAACT 91) 
Jmeet + 1 
G@eTO 103 
moe DOPlift |! = 196 
i= Of 
ml GGt(!l) = 0. 
CCt1) = GI1)/AA(1 91) 
DOwiwsS 6€6fld6= Ba 
LIM3 = ait | 
Wememre K = Lok ims 
eee oer) = CCCT) FF CAACI SK) * CCCUK) ) 
Goal) = C( 1) =aCc (7 
Meo CC(L) = CCCI)FBAUTs1) 
X(6) = CC(6) 
Demis it = ise 
Il = 6 I 
Chee = IT + 7 
DO 114 K = LIM496 
114 X(II1) = XCTT) + (AACTIsK) * XCK)) 
mee ACIT) = CO@MTyY = xXt(IT) 








GO TO 118 
PRINT 117 


PORMAT (2 8)H SINGULARITY —- NO SOLUTION) 


RETURN 

END 
Seer ouUTme OR STORAGE OF THREE SETS OF VALUES 
OF AXTAL '.OAD AND THE FOUR DEFLECTION COEFFICIENTS. 

LABEL 

Ep esa is: 

FORTRAN 

GWBROUTINE STORE (AXs Ws Zs Ss LL) 

DERMENGION W(6)s 2(3)5s S(306) 

ZC) = AX 

DEMOS =ahs6 

Suess i) = WT) 

RETURN 

END 
Pans SUBROUTINE IS JUST W SHIFTING PROCEDURE 
OF VALUES FROM ONE STORAGE LOCATION TO ANOTHER. 

LABEL 

Hiccu ies 

FORTRAN 

SUBROUTINE BAGH (Zs S) 

BIMENSTON Z2tay75 S (356) 

DO 2904 I = 196 


DO™ 204 J = "Faw 
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SC Ie!) =aeSakee Agel ) 
Pee206 Iieegloe 
aS = 1) + F 
Orit) = Z2(J69 
DO 206 JI = 196 
mrosJl) = 6. 
Z(3) = O-6 
RETURN 
END 
SUBROUTINE FOR PRINTING 
LABEL 
LIST8 
FORTRAN 
SUBROUTINE CPRNT (PNOR'Ms AXs Ws ERROR» RONE) 
DIMENSION W(6) 
SUM=W( 1) #*#24+W(2 )#X¥24W(5) F¥249 0% (W(3) HK D4W(G) HZ) 4+25ekW(G6) H*D 
SUTRAINT= AX 4°°341416%*2 *°SUM/(8. # RONE) 
PRINT 2005 PNORMs AXs (WfT)s IT=196)9 STRAIN» ERROR 
MOnwnl (PS.29F11 05> 6F12665 F115» F127) 
RETURN 
END 
SUBROUTINE FOR COMPUTING ALL THE NECESSARY 
CONSTANTS IN THE MAIN PROGRAM 


LABEL 





Pos 
. FORTRAN 


SUBROUTINE CONST (Us ALPHAs PNORMs Rs Bls B2s B3) 


S 
ei eve 


ola © wo 
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DIMENSION P(13)s R(54)5 B1(4)s B2(5), B3(3) 


Par) AL PHA#*2 


iN 


RUZ) = 1.7 PUES 

Peo o™ 1 = ae7 

eiSe (1 -— 2)8e?7 

mie PSP (T)/ 0PT1) + Y1)**2 
DOwP301 J = 8511 

yae= (oY) - 67EFZ 


Se = Pewee + YA * Pll] ) }e2 





P(i2} = P(1y7t4. + 94) ¥#2 
P(13) = P(1)/(4. * P(1))*#2 
R(T) @= Pr? 

R(2) = P(2) 

R(3) = P(3) 

Ree Se P(e) 

R(5) = P(8) 

Rae = P(1) P(2) 
R(7) = P(1) 81.*P (2) 
RG = P(l) + 625.*P(2) 
R(9) = 81.6#P(1) P(2) 
R(10) = 625e*P(1) P(2) 
Rvil) = Pig: ) P(8) 
R(12) = Pid ) 81.*P (8) 
R(13) = P(2) P(4) 
R(14) = P(2) 81.*P (4) 
R(1%) = P(1) 3.%*P (4) 
R(146) = P(2) 34%P(8) 
R(17) = 3e*P(1) 64, %*P (3) 
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R(18) = S4RP2) + 6eeFPES) 

R(*9) = G6@e#P(°) + 2746, *P 64) 

R(@O) = GeerP(3) + 274.*PT8) 

R Gas) = P(4) + em ) 

R(22) = Prt.8.)) + oS) 

R(23) = PCa) & P(9) 

R(24) = Sst Poe] pet CS FP(3) *% 25e*P(8) 

R(25) = PCePR C2) + 64.e*P(3) + 25e*P(4) 

R(26) = turret + 44] ¢*P(8) + 256e*P(9) 

R(27) = Comae G2) + 44] .%*P(4) b 256e*P(5) 

m2) 2 fame P C3) + G2ZoeeP (a4) + 6250" P (8) 

R(29) = 9e*P(1) + 646*P(3) + 12256*P(4) 

R(30) = 9e*P(2) + OHeFPO3) + 12252¢*P(8) 

Rvs) = 25erRta) + 64.2%*P(3) + 92*P(10) 

R(32) = 25% (2) + 64.*P(3) + 9e*P(6) 

BASS )>> 64e*P(3) + 44] .#P(4) + 169e¢*P(13) 

R(34) = 64e*P(3) + 441] ¢*P(8) + 169e¢*P(12) 

R(35) = 2omteee ) 6+ Seem 57 t)6| 6BAB Ie EP (13) 

R(36) = oO ) + 516e#P(3) + 8281-*P(12) 

R437) = 155200) + BOSSOMEPIT2) + 2856le*#P 013) 

R(38) = 9h) = 4e*¥P(2) + 166*P(3) + P(8) 
R(39) = 4e*P (8) + Ge*P(2) + Ge ¥e yo P(4) 
meee = 366*PI1) + S6eRWP (2) + 2256%P(4) + 22506%P(8) 
Rel) = 25emP C1) + 4e*P(2) + B8le*P(8) + 1d5e6*P(9) 
R(42) = Liar ae) + eer (2) + B1le*#P (4) + 166*P(5) 
R(43) = 36e¢*P(1) + Pare je ae?) 9e*P(13) 
R(44) = PS i 36e*¥P(2) + 144e¢*#P(5) + O@FP( 12) 
Res) “= cemenl) + 16e«*P(3) + 4Ge%#P(4) + 92*P(8) 






if 

| | Le ° 
—e ; Peter bas - | st omen e +i Gare 
ie 


shar, t — | 









y aif ‘ower wT 














i) tear 
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R(46) = 2ePR(2) + herr (3) + 9e*#P(4) + 49e%*P(8) 

RG47) = 16eP13) + Ble#P(4) + P(8) = 8] eer (10) 

RGSS) =" eeerrrs) + Eee) + Came Pr (6) + B1e*P(8) 

R(49) = emery +t 22562) +t 256 etP (3) + PC6) 
Re) oe er cl) + eee) + 256e%P(3) + mae ly Oi, 

pee a Gee (3) + 2401 6 FP4) + 49096-*P(5) + 2401-*P(12) 

pemere 256. FP (3) + 2401 oP (8) + 4096-*P (9) + 2401.#P(13) 

R(53) = 2.*P¢1) + 648*#P(3) + 2566*P(5) + 2256*#P(8) + 49-#P(12) 
R(54) = 2e*P(2) + 646*¥P(3) +) 22506*P(4) + 2566*#P(9) + 496%P(13) 


DO 302 i898 @¢43 
m3 = 2 °° te 1 
Bl(l) = 2564 * P(]) **PNORM/(Y3 * 3461416%~«6) 


= iar ee 1 Beh) /(3.1416%*2) 


oon | ) 


B1l(4) 


B27 

po 303 J =4F95 

ol 2 ee) Pe 

B2( J) =6e07 Uses (ls UF 2 ee P11) ) 

RETURN 

END 
SUBROUTINE FOR COMPUTING THE COEFFICIENTS 
FOR THE USE OF THE CROUT REDOX METHOD 

LABEL 

L1STs 

FORTRAN 

SUBROUTINE SESAME 

COMMON Ws AXs Bls B2s B3s Rs MGs As C 

DIMENSION W(6)s 81(4)s B2(5)s B3(3)5 2154) Al6s6)5 Cl6) 


GO TO (491+ 400}+ MG 
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moO Alls1)=B3 01 pew?) 
A(291)=B3(1)*W(2) 
A(3s1)=B3(2)*W(3) 
A(491)=B3(2)#W(4) 
A(591)=B3(1)#W(5) 
A(691)=B3(3)#W(6) 
GO TO 402 


401 Ali191)=82(1)4B3(1)*AX+32.*R(6) *¥W( 1) **¥9-6.*R(2) #W(1) *W(2) 


1 —6e*R(1)*W(1) W132) 4R038) HW 2) #*#24+R639) HW(9) HD 

Zz +32¢MR(3)#W(2)#W(3)+81.*R(21)#W(4)K*2-18.6*R(11)#W(2) WG) 
3 —18-6*R(13)*#W(3)#W04)-66*R(16)#W(2)*¥W(5)-66*R( 15) ¥W( 3) ¥W (6) 
u +162e*R(5)#W(4) W155 )4+1620*R(4) WIG) #W(6)4+R( 41) #W(5) *¥2 

5 +R(42)*W(6)**2 


AO2s 1 )= se eR OR Oe Dees 72 errR( 38) AW) HWI2)4+32-6*R03) WO 1 PW) 


i —18e*R(11)*W01) ¥W(14)-R( 25) *¥W( 2) #*2-2,.*R( 24) HW 2) FW 2) 
? +R(40)¥*W(3)*W(4)-64*R(16) #W01) #WI5 74+ 26*RI18) HW 2) HW 5) 
3 +R (20) *#W(3)%#W(5)4R(19) #¥W(3) ¥W(6)-144.¥R( 23) #W(4) FN(5) 
4 —-9.#R(45)#W(4) *#W(6)-R( 27) #W(6) *¥2 


ee yas. men (CP) eee og eR Cay eW01) eWI2) +206 RI 39) FW 1 ee 8) 


] [198 FER) HWE] ) FW A PERODS PW 2D PH 2—2.ER( 25) %W( 29 HW B ) 
2 +R(40)#W(2)*W(4)-66*R(15)*W(1)*W(6)4+R( 20) *W(2) #W(5) 

3 +R(1L9O)*#W(2)HW(6)F20*#R( 17) EWL3) HWI6) -FOe*®RI G46) XWI4) FW 5) 
4 —-144.#R(22)*W(4)*#W(6)=R(26)¥W( 5) #¥*2 


AC491)=-1238¢MR(11 )#W(1)*WO2) -184*R (13) #W01) eW03) 


] +162e#R(PT)*HW(1)#W(4)4R( 40) ¥W(2)*W( 32) 4+16208#R(5) (1) FW(5) 
2 +162e#R(4)HW(1)*W(6)-144e*#R( 72) *W( 2) EWS) 
5 —9.#R(45)¥#W(2) #W(6)-92*¥ R146) FW( 3) HWI5)-14464*R( 22) FW(3) FW (6) 


4 +18 e*R(14)*#W(4) #WI5)4+182*R(12) #W04) WG) 





7 
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. 
= 7 - = - 
> 7 ‘ ; ~ 
j ‘. - ~ eo 
= 4 ue 
* | § _ j - rm 7 
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= pee eiiee 

- opecvsiiiat P mate 
| reir” 
oie wowelt =] rt jet sprt- 


Fe views ial (Pye oh mart tee ", reine 
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5 +729 e#RUCLIREW(5) EWE) 


A(591)=~6e#R(16)#W(1)*#W(2) 4167 o*RI5) W101) FW(4)+R(18) FW(7) H*D 


] +R(20)*#W(2)#W(3)4+90*R(14) *W(4) *#7-144.*#R(23)*W(2) HW (4) 
2 —9.%*R(46) *W(3) #W(4)420*R(141 ) W001) *™(5)-20%#2( 26) N13) *W(5) 
2S +729 e¢*R(21)*#W(4)*#W(6) 


A(691)=-6e*#*R(15)*W(1) #903) 4+162.*RV4)4N (1 HW(4)4+R(17) #W( 2) # #2 


1 +R(19)#W(2) W139) 492 #R( 12) 4W(4) F¥2—4 Ce *RI(G45)HW(I D2) WG) 
2 —144-6#R(22)#W(3)*#W(4)42,*R( 42) ¥W(1) HW(6)—-22*RI27) HW(2) N16) 
3 +729e#R(21)*#W(4)#W( 5) 


402 Al1s2)=-3e*R(2) *W( 1) ¥¥#242.*R( 38) *WE1) 4W(2)4+326*R( 3) *W( 1) HW 3) 


—-18e¢*R(11)*#W(1) #W04)-2.*R124)#W(2) *W(3)-R1 25) HW 3) # HD 
2 +R(40)#W(3)*W(4)-66*R( 16) *¥W(1)*¥W(5)422*R(18) *W(2)%*W(5) 
3 +R(20)#W(3)*¥W(5)4R(19)*#W03) ¥W(6)-1446*R( 23) WIG) HW(5) 
& —9e#R(45 )*#W(4) #W(6) -RI27) #W(6) **2 


A(2s2)=B2(4)+B3(1)*AX+R( 38) #W(1)**2-26#R( 24) *#W(1) #W(3) 


i +3e*R(9)#W(2)**¥24R( 28) ¥W(3)*¥*2-4GB66*R( 1) FW(2) HW 4) 
2 +R(43)#WUA) HH 24+2.*R1 18) WOT) AW05)-5120FR09) HW(3) MWD) 
3 —2e*R( 29 )¥W(3) #W(6)4R050)4W(5) #¥24+R(51)4W(6) HQ 


A0392)=160*R( 4) OW (1) FF 2—2 0 FRU 2G) FW) HW 2) - 20 FRO25 EW AW  ) 


1 +R(40)#Wl1)4#W(4)42.%R( 9289) 8W(7)4W(2)4R(720)%#N(1) #W(5) 

2 +R(19)#¥W(1)*¥W(6)—-5126*#R13) *WI7) #WI5) -2e*®R( 29) #W(2) ¥N(6) 

3 —2e*R(30)*#W(3)#W(15)-5122*#R(13)%W(9) 4 (6)4+92#R(54) W104) FW (5) 
4 +Ge%#R(53)#W(4)¥N (6) 4656 0*R(3)#W(5) #W(6) 


Al4s2)=-9e*R(11)#W(1)*#*24+R (40) #w(1) W093) -2 492 FRI1) FW 2) HQ 


it +1 8s Rimes Rew? peew04 )-144.*R( 292) *#W01) twW05 )-9 oF RUB] eS) eee 
2 —9e*R(45 )*#W(1)#W(16)490#R(54) W139) HW(5)490*FR(53) ¥N(3) RWG) 
B) —9e*R(33)*W(5) ¥W(6) 


A(5927)==30R(16 PEW] 1 F942 HR (ADH VC EEC 2) FRO) ee ea 





| yen en nero a FIRS 
» os 
4 


2 
a. (=o : 
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a § 
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. Cd | 7 
| pecan Sy eee 
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‘@@,44 
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Al69?)= 
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Aves 3) = 


Mi4s3)s 


Pa 593 ) = 
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=1446#R( 23) ¥W(1)*W(4)-512.4#R(2)HWI2) HW(3)—-P (139) (3) HED 
+#9e%#R(54)#W(3)4#W(4)426%R(50) ¥W(2) ¥W(5) 46560 *¥ R13) HN (3) *W(6) 
—-18e*R(313*#W(4) *W05 )-96¥R(33)#W(4)4W(6)-R( 36) *W(6) ¥*2 
R(19)*W0] )#W03) ge FRIG45) W001) WIG) -2 RI D9) H¥W(D) HW) 
—-2560e#R(3)#¥W( 3) #249 24R (53) 4W(3)4WI4)-2.*#R( 27) #W(1) W156) 
#2 eo#R(51 )¥W(2) #W(6)4+6560*R13) W139) FW(5)-92#R( 33) HN(G) HW(5) 
—2e#*R(36)*W(5)*¥W(6) 

—-3.#R(1) ¥W01) ¥*#2422.#R( 39) WI] ) HW(3)43706%*R(3)4*W(1) *W12) 
-186#*R(13)*W(1) #W04)-RU24) 4W(D) H¥D—2.¥RI 25) HW(2) HW 3) 
+R(40)*¥W(2)¥W(4)-65.*R(15)*¥W(1) #W161420*R(17)4¥N(13) FW(6) 
+R(20)#W(2)#W(5)4+R019)#W02 )*W(6)-Fe#*R(46)HW(G) HW(5) | 
=144.#R(22)#WI4S¥W(6 )—-R(26) *W(5) #¥2 
16e¢*R(3)#W(1)¥*¥*¥2-2,%R( 25) WIL) HN (37-2 FRI 24) KW 1) FWD) 

+R (40)4#W(1)#¥W(4)42,.#R( 28) %W(2) #W(32)4R(20)4#W(1) #W05) 

#2019) *#W(71)#W(6)-2,.*R( 30) FWA) KW(5) -51224*#RI3) WI 2) 4W(4) 

—5 12 e#R(3)HW( 2) HW05 )—2eFRI ZI) HWP? )D*FW(6)F4+90#R( 54) ¥W(4)HN(5) 
+9e#R(53)#W(4)#W16)4+6564*R( 3) #WI5)#W(6) 
B2(2)4+B3(2)*AX4+R(39) W011) #H 2-7, #R( 75) HWI1) HW02) 
+R(28)#W(2)**¥24+3.4*R( 7) W038) *¥¥O—-GB6 RI 2) FW (3) ¥W(G) 
+9e*R(G44)#W(4) #*¥242,.*R117) FW(1) X¥W(6) -2 2 FR 30) HW( 2) FW(5) 
—-5126*#R(3)#W(2)*W06)4R052) ¥W(5) #¥24+R(49) HW(6) #¥2 
—-9-.#R(13)*W(1)#*¥24R(G0) #W(1)¥*¥W(2)-24326%#R(2) ¥W(3) **2 

+18 eMR(44)¥W(3) #W14)-9-6*#R1 46) ¥W(1 W105) -9 FRI 32) HN (6) HD 
—-1446#R(22)*#W(1)#¥W(6)4+9.*#R( 54) HW(2) HW(5)4+92.*R(53) FWD) HWS) 
-9.*¥R(34)¥W(5) ¥W( 6) 

R(2O)*#W(1YHW( PD) HDeHR( 46) WIT) HWIG)-25EA CERI 3) HNI 2) RED 

=2 eo FRI 30) #W(2) #W(3)496*R(1 54) FWD EN 04) HDeFRI PG) HW] YEW ) 


+PeBER(SPIEW(ALEW(5) F656 FRI 2ALHWI AD HN(H) HOP FRI 324) FH(C4) HM (6) 









a ea 


erm . 
ri e°oeret : : J hee yom re ere ah. |? 
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A( 1 o4 )= 
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—2e*R(35)#W(5)#W(6) 
—3e.*#R(15)¥WI1)*¥¥*¥242.*#R(17)¥W(01 413) +R119) HW) NZ) 
—1446#R(22)*W(1)*W(4)—-R(29) ¥W( 2) **#2-5122*R(3) Hw (2) HW) 
+9e#R(53)#W(2)#W(4)4+20#R(49) WI 3) 4W(6)4+6566.*R( 3) *N(2)%W(5) 
=Our mn (34 )eW(H) FWD) -18.*RI3B2)4#W(4)eW(6)-R( 35) ¥W(5) #2 
162.MR(21)*W(1)*W04)-18.*R(11)*W(1) *W( 2) 
—-182*R(13)%*W(1) ¥W(3)4R( 40) FWD) ¥W(3) 4162 064R(5) *W(1) HN(5) 
+1620#R(4)*W( 1) #W16)-1446*R( 29) W022) *#W (5) 
—9-#R(45)*#W(2)#W(6)-9-¥RI46) FW 3) KE (05) -144.*#R122) HW (3) H#N(6) 
+18 6¢*R(14)%*#W(4)*W(5)4+182*R(12) ¥W(4) ¥W( 6) 


—729e¢#*R(21)*W(5) #W(6) 


A(294)=-9e*R(11)*W(1) **#24+R (40) #W(1) FW 3)-243.*R(1) *W( 2) **2 


+20#®R(43}*#W(2)*#W(4)-144,*R(23)4W01)XW(5)-9eFR(45) WI 1) RW 6) 
+9e#R(54)#W(3)*W(5)4+90*R153) WI 3) *W(6)-92*#RI 31) HWI5) H¥2 


mo ome USS EWS yEWS) 


A132 s4)e=9.*#R (139 W (1) e422 +R OD EW 1 EW? )—- 9436 FR(02) ¥W13)4*2 


Pol» 4 ) 


— 
= 


+18 e#R(44)#W(3)*WI4)-9.*RI4G6)HWI1 DW (5)-96*R( 34) HW(5 ) HW(6) 
=144.*#R(22)#W( 7 SWIG) F9,#RI54) MWD) HWI5)4+92*R( 53) HN(2) RWG) 
#9. *RCUS2) RWG ) * FO 
81e*B2(1)4B3(2)*AXA4+81 o#R( 21) W111) F249 0*#R( 43) HW( 2) HH 

#243 e#R(6)#W(4) ¥¥24+967*R( 44) €W( 3) FK2418 o*R(14) WI 1) ¥W(5) 
+18eMR(12)*#W(01)*W16)481e*R(47) ¥W(5) #2481 4%#R( 48) HW(6) ¥¥D 
+2592.*R(3)*W(5)#W(6) 

B1e*#R(5 )*#W(1)¥*¥2418.*¥R114) FW) HW(6 49 e*R( 54) FW?) EW) 
—144e¢*¥R(23)#W(1)*#W(2)-9.*R( 46) ¥W(1) *W(3) 

+729 e*R(21 HWI 1) HW (6)—-1B 6 *¥RI351) *W(2) *#W(5) 
—9e*R(33)*W(2)*W(6)-9.*#R1 34) HW( 3) EM ( 6) +162 0% R147) FW 4) N05) 


+2592 6#R(3B)#W(4)#W(6) 
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A(634)=817.#R(4) *W(1)*#*2419.*R(17) W071) H¥W(4) -9-*R( 45) AW) HD ) 


1 —-14446#R(22)*W(1)#W(3)4+9.*R(53) HW 2) HW 3) 

E, +729 e*R(21) ¥W(1)#W(5 )—-9 eo *¥R133) HWI2) WI 5) -9e*R( 34) HW(3)%W(5) 
3 —186%*R(32)#W(3)#W16) 42592 6# R12) HWI4) HWS) 

+ +1626#R(48)*W(4)*W(6) 


N95 SG eR (PE VFWIT )AW02)4162-*R(5) HW) W004) + 2eHR (4) FW (7) Fw 5) 


1 +R(18)*¥W(2) #8 24R(20) #W(2) ¥W13)-144.*R( 23) %*W(2)%*W(4) 
2 —9.%*R(46)*#W(3)*#W(4)4+9.*R(14)#W(4) #¥ 2-2 2*R( 26) %¥N(3) HWS) 
3 +729 e#*R(21L)#W(4)#W(6) 


A(295)=-3e#R(16 )¥W( 1) #242 -*R(18) FW(1) HW 2) 4R(0 20) *W0 1) HW 3) 


] -144.#R(23)*W(1)*#W(4)-R( 30) #W( 3) #¥2-5126*#R(1 3) HW(2) HW (3) 
? +9e%*R(54)#W(3)#W(4)42.*R150) *W(2)*W(5)+6560*R( 3) *W(3) eW(6) 
3 —186*R(31)*#W(4) #W15)-96#R( 33) *W(4)*W(6)-R( 36) WIG) **2 


A(395)=R(20)#W(1)#W(2)-9e#¥RI46) #W(1)#W(04)-2.*RI 26) WI) ¥W5) 


1 —256e#R(3)#W(2) ¥¥#2-26*#R( 30) AW 2) W193) 4+90*R(54) FW (2) RWG) 
2 +20%*R(52)*W(3)#W(5)4+6566*R(3)#W(2)*¥W(6)-9e*¥RI3B4) EW(4)*W(6) 
5 —2e*R(35)*#W(5) *#W(6) 


A(425)=81-e*R(5) *W(] ) #*¥#2=-144.%R(23)*W( 1) HW(7)-9-¥RU 46) FW (1) FW ) 


i +18 6*R( 14.) *¥W(01)#W04)4+7292*#R( 21) Wil) HW(6) 

2 £9 0#R(54)#W(2)*W(3)—180¥R(31)¥#W(2)*W (5-94 ¥*R( 33) 4002) #W(6) 
3 —9e*R(34)*W(3)*¥W(6)4162-e*R(47) WIG) FW 5) 

ty +25922#R(3)*W(4) #W(6) 


A(595)=B2(5)4+B3(1)*AX+R(41) #W(01)*#2-2,*R(26) #W(1) #W03 ) 


1 +R(50)*¥W(2) ¥*¥242(52)*#W(39)**2-178.*R( 31) HW 2) *W(4) 
2 +81e*®R( 47) *W( G4) #¥#2-2 .*R( 35) FW( 3) WIG) +30 FRI 10) HW(5) HHO 
3 +R(27)#W(K)*#? 


A(695)=729e*R(21 )*W(1)*#W14)4656.*R( 3) ¥WI2) HW 3) 
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